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ABSTRACT 


The purpose of this research is to develop algorithms for mapping terrain 
characteristics from a 100 meter grid representation to arcs and nodes of a 
transportation network. The algorithms capture elevation and trafficability parameters 
as they relate to both the arc itself and to the off-arc characteristics. These network 
parameters are directly usable in appropriate optimization algorithms to determine 
minimum travel time path and minimum distance path through the network for a 


variety of maneuver unit types and missions. 
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I. INTRODUCTION 


A. PURPOSE AND GOALS 

The purpose of this research is to develop algorithms for mapping terrain 
characteristics from a 100 meter grid representation to arcs and nodes of a 
transportation network. The aigorithm must capture elevation and trafficability 
parameters as they relate to both the arc itself and to the off-arc characteristics. These 
network parameters must be directiy usable in appropriate optimization algorithms to 
determine minimum travel time path and minimum distance path through the network 
for a variety of maneuver unit types and missions. 

After. reviewing the literature and discussing the problem with several terrain 
model experts, it was determined that no automated procedure for accomplishing this 
mapping exists. Therefore the procedures and associated algorithms presented in this 


thesis represent original research and development. 


B. AIR LAND RESEARCH MODEL METHODOLOGIES 

The purpose of the Air Land Advanced Research Model ( ALARM ) is to 
develop a systemic corps level combat model, that is, one that operates without human 
intervention. This idea is not new, but ALARM takes a different approach than 
previous efforts to simulate the military decision process. 

“Іп most combat models, the structures for planning and for execution are 
included in the same program. In ALARM the planning function and the execution 
function will be separate programs. The Planning model will be a carefully designed 
collection of programs that provide orders to the Execution model at the appropriate 
times. 

A second major difference between the ALARM design and existing models is in 
the method used to make decisions. All of the models reviewed when developing the 
ALARM strategy depended either on threshold parameters or on a series of “IF 
THEN” decision rules to make decisions. 

Combat models using a threshold parameter strategy usually require a large 
number of threshold parameters or groups of parameters. In such models a threshold 
parameter ( or group ) is needed for every possible decision that any unit represented in 


the combat simulation might have to make. This leads to numerous parameters and 


thresholds. Most parameters of this type have no reasonable real world counterparts. 
In other words, parameters used for thresholding are frequently contrived and based on 
unexplainable and undefendable measures. This makes it difficult to create data bases 
for such models, and difficult to place any confidence in the analytical results from the 
models. 

Models that use a series of decision rules are lacking in that they require a very 
limited view of the decision process. A decision can be made only if a series of decision 
rules to support the decision exist within the software that is developed for the model, 
and if those specific rules are satisfied. When a situation is encountered for which there 
15 no specific decision rule, a default rule must be implemented. These default rules are 
frequently inappropriate. The more complicated the decision rules become, the better 
the model usually is, but the more difficult it becomes to maintain the model and the 
harder it is to analyze results from the model. The cause and effect relationships 
established when running the model are usually closely tied to the decision rules. Thus 
the output or results from such a simulation are built into the model. The foregoing 
discussion emphasizes problems with the methods currently accepted and used to 
model the decision process. These methods are recognized because they have their 
uses, applications and foundations in the human decision making process. Human 
beings do, on some level, use the threshold initiated decision rule process in their 
. decision making procedures. These concepts are usually used to start or constrain the 
decision process and not to make the decision itself. 

The premise of ALARM is to use both threshold strategy and decision rules, but 
not within the same framework as other models to date. Thresholds will be used to 
determine when the planning or decision making procedures should be executed. 
Decision rules will be used to limit alternatives. They will not be used explicitly to 
make combat decisions. ALARM will use network based methodologies to itemize 
alternatives so that a decision can be made. It is believed that this methodology 1s 
closer to what actually happens in the human decision making process. 

ALARM decision processes, called Decision Tasks, will use three unique 
methodologies that are the mainstay of the ALARM concept and make it different 
than other current research. These methodologies are time domain networks, cartesian 
space networks, and the Generalized Value System. Each methodology is briefly 
described below. 


1. Time Domain Networks 

Time domain networks are designed to handle the planning function or 
activity within the ALARM system. А time domain network consists of nodes and arcs 
( or links ) connecting the nodes. The key to time domain networks is that arcs do not 
represent distance, but represent the passage of time or the completion of a sub- 
activity which leads to the completion of the entire activity represented by the network. 

The time domain network is used to develop high level mission requirements 
for all subordinate units. It is similar to a PERT or CPM network in that a sub-activity 
cannot be started until all merging activities are completed. The input to the activity 
network is a desired actıvıty completion time and an acceptable friendly force attrition 
level, and may include a required level of attrition to enemy forces. 

2. Cartesian Space Networks 

The Cartesian space network is a series of networks each representing a 
different aspect of the battlefield. Each network represents physical connections 
between points on the battlefield. ALARM will have three or more Cartesian space 
networks. The three networks that have been identified are terrain and transportation 
networks, communication networks, and logistics resupply networks. This thesis deals 
explicitly with the terrain and transportation networks. 

3. Generalized Value System 

The generalized value System (GVS) is a £i for quantifying the 
capabilities and importance of entities on the battlefield at some future time (t-- x). It 
does this by developing algorithms to predict future entity or situation states at time 
(t+x) based on the situation at time (t). Thus GVS provides the framework for 
forecasting future states of entities in continuous time. 

Various combinations of exponential functions are used to represent these 
forecasts. For example, the power and value of non-combat entities such as bridges 
and road interdiction points; combat support units; service support units; and combat 
units can all be determined from a continuous function with respect to time using the 
various exponential algorithms. These values can then be used to determine the specific 
place or target against which to plan an air strike or a ground attack at time (1+ х). 
Thus different targets can be selected by the Decision Task depending on the time at 
which strike assets are available. The decision will be based not on the current value of 
the target, but on its predicted future value in relation to the rest of the Air Land 


battle. To do this, all targets and units must be assigned consistent values so that 
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comparisons can be made across the diverse set of targets available on the Air Land 
battlefield. The common bond among these varying targets is a value system that 
measures the target’s value in terms of its usefulness to combatants. Thus the value of 
a bridge is not measured simply by its width or length, but is tied to the number of 
enemy or friendly units that could cross the bridge to maneuver into favorable combat 
positions. As time passes, the bridge’s value could increase based on heavy combat 
areas and then could decrease when all relevant units have successfully crossed the 


bridge. 


C. SCOPE OF THE THESIS 

| In attempting to achieve the goals of this thesis, various terrain modelling 
techniques will be described in Chapter II as a motivation for the network algorithms 
developed. The algorithms for a single arc will form the basis of the total model and 
will be presented in Chapter III. The fuil network implementation wiil be presented in 
Chapter IV to determine the minimum time and minimum distance paths on an 
undirected graph. Chapter V will discuss the results of several runs of the model to 
illustrate its various capabilities. Finally, Chapter VI will provide recommendations for 


continued enhancements of the model. 
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Il. BACKGROUND 


A. TERRAIN MODELLING TECHNIQUES 

Representing combat functions as a system of networks is a departure from the 
current practice of explicit terrain representation. Current modelling techniques require 
large data bases to represent numerous aspects of terrain such: as slope,. elevation, 
forestation, and cities. The network approach is an attempt to. represent. the terrain, 
primarily through use of the transportation system, as an abstraction of its actual state. 
The topology characteristics required by explicit terrain. models will be: reflected in the 
arc and node attributes. 

The research. model will gain an economy of representation of the topology 
because it has the flexibilitv to represent only those terrain features applicable to a 
stated level of resolution. The corps rear is sparsely populated with combat and 
support units relative to the overall size of the corps sector. Furthermore; large 
portions of the corps sector will never support unit movement or facilities and have no 
reason to be modelled. A network representation allows one to. select’ only those 
features which, because of their nature, give the controlling force a marked advantage 
over the other force ( e.g. key terrain ). Additionally, those roads or cross-country arcs 
which can either expedite or delay movement can be identified and modelled. 

In a corps level model, direct fire combat requires observation, detection, and 
lines of sight along a narrow strip of terrain relative to the overall size of the sectors of 
the two forces. Existing methods of terrain representation may still be required to 
model the details of such combat. The majority of units moving in a corps sector, 
especially the CS/CSS assets, do so without receiving direct fire. The movement of the 
units can be modelled via a network abstraction of the transportation system. 
Numerous algorithms exist which allow one to exploit the features and. structure: of 
such a mathematical model. 

Many of the terrain modelling techniques in current use restrict the model to one 
level of resolution and force the terrain data base to store attributes for large sections 
of the battlefield which are not used. To put current terrain modelling procedures in 
perspective, the following is a short discussion of three terrain modelling methods: hex 
terrain ( used in CORDIVEM ), digital terrain (DYNTACS) and network model 
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В. HEX TERRAIN 

The Corps-Division Evaluation Model (CORDIVEM) is a corps level model in 
which the terrain is represented as a series of hexagons, (hexes), each approximately 3.5 
kilometers in diameter. The terrain underlying each hex is characterized by a three digit 
array which represents urbanization, forestation, and roughness. Each individual 
attribute is an average of the Defence Mapping Agency’s data base of points 12.5 
meters apart which underlie each CORDIVEM hex. The attribute based on the 
resulting average is then assigned a code in CORDIVEM and assumed homogeneous 
throughout the hex [Ref. 1: p. 16]. Figure 2.1 illustrates the codes for each hex 
attribute. 


ROUGHNESS OADS 

0 - 15% 0 - none 

- 304 l = tertiary 
- 504 : 2 - secondary 
- 504 | 3 — primary 


1 
2 
3 


RIVERS 
0 - none 
1 - stream , 
: | dec river URBANIZATION ———1 
О = 154 
302 
502 


502 





Figure 2.1 CORDIVEM Terrain Features. 


The geometry of the hex is used to represent unit movement, roads, rivers and 
obstacles. Unit movement and a road system, if one exists, are modelled from the 
center of one hex to the center of an adjacent hex. Movement on a cross-country route 
is penalized by reducing a unit's allowable speed. Rivers are modelled along the edge of 
a hex, and coded as indicated in Figure 2.1. An implicit bridge is assumed at the 
intersection of all roads and rivers. The bridge can be explicitly modelled if it may 


affect movement factors or possibly be destroyed. An obstacle occurs on the edge of 
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the hex and is explicitly inserted into or removed from the model. All obstacles are 
created before the execution of the model; the dynamic construction and reduction of 


obstacles are not modelled. 


C. DIGITIZED TERRAIN 

The Dynamic Tactical Simulation (DYNTACS) model is a two-sided, dynamic 
Monte Carlo simulation capable of modelling the combat process from the individual 
crew to battalion engagements which utilizes a concept referred to as digitized terrain 
[Ref. 2: p. 9]. Macro terrain in DYNTACS is usually represented as 100 meter squares 
with attributes for the corners being supplied by data available from the Waterways 
Experiment Station. These corner attributes provide information about such features as 
elevation, forestation, and location. Each square is then divided diagonally, resulting in 
a series of equal sized triangles, varying in slope and orientation. 

This approach assumes that the terrain modelled by each triangle 15 
homogeneous and that elevation changes spacings. Sudden changes in the terrain or its 


surface are not detected; instead, they are represented with geometric overlays. 


D. NETWORK MODEL FOR MOBILITY IN ALARM 

The Cartesian space network compacts the terrain database still further by 
limiting battlefield movement to a subset of battlefield locations. The allowable 
locations are represented as a network overlayed on the battlefield map. Nodes of the 
network represent physical locations on the map. Arcs of the network represent 
possible movement paths between nodes such as roads, trails, or likely cross country 
routes. 

Each node and arc of the network has mobility attributes which determine the 
limiting speed of combatants moving through that node or along that arc. The 
attributes which are stored in the network representation can be either terrain features 
( which will support an explicit mobility model ) or mobility multipliers ( for a 
simplified implicit model ). 

Added flexibility can be gained by storing several sets of mobility attributes for 
some arcs. The main attribute set might represent movement along a road, while a 
second attribute set could represent cross country mobility in the near vicinity of the 
road if the road itself cannot be used. 

The main problem with network models for battlefield movement is that 


movement is restricted to the network, and this may limit tactical flexibility. 
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Preparation of the network must be carefully coordinated with the tactical scenario to 
ensure sufficient movement opportunities. This limitation would seem to be most 
severe in the direct fire battle zone where cross country maneuver is common. Areas 
away from the front lines can more readily be represented by a network model. A 
major advantage of the network model is that it allows efficient network optimization 
algorithms to be appiied to the problem of route selection for moving combatants. 

To date, data for network attributes have been derived manually using map 
analysis. This technique is both extremely manpower intensive and inaccurate. The 
requirement for an automated process to generate network attributes directly from 
digitized terrain data bases serves as the motivation for this research. The algorithms to 


determine attributes for a single arc are described in the next chapter. 
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III. SINGLE ARC METHODOLOGY 


A. OBJECTIVE AND PROBLEM DEFINITION 

The objective of this research is to develop algorithms for mapping terrain 
characteristics from a 100. meter grid representation to arcs and nodes of a 
transportation network. The algorithm must capture elevation and. trafficability 
parameters as they relate to both the arc itseif and to the off-arc characteristics. These 
network parameters must be directly usable in appropriate optimization: algorithms to 
determine minimum travel time path and minimum distance path through the network 
for a variety of maneuver unit types and missions. 

Methodology for the single arc attributes which is used for the Cartesian space 
network is presented in this chapter and is related to the computer program for the 


single arc attributes in Appendix C. 


B. MODEL DATA FILE 
The program used four disk data files, one for 100 meter grid square data, one 
for the node characteristics, one’ for the arc characteristics which represents the 
networks overlayed on the gridded terrain, and one for the speed data related to 
Maneuvering units. 
1. 100 Meter Grid Square Data 
The grid square data represents the altitude and characteristics of each point 
spaced 100 meters apart on the terrain. The source of the data is the Combined Arms 
Center at Fort Leavenworth, which is the same data base originally used for : 
CORDIVEM. The gridded terrain is used to determine the attributes of networks 
overlayed on the terrain. A typical network overlayed on the digitized terrain 1s shown 
imn Figure 31 
Two characteristics are used to describe a point of gridded terrain: three digit 
altitude and an integer surface feature type. The altitude is the height from sea level or 
simple map altitude. The characteristics of the surface feature type in the network are 
used for determining the arc width for maneuver along the arc. Table 1 shows how 


data for the 100 meter resolution surface feature types were encoded. 
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ТАВГЕ 1 
CHARACTERISTICS OF SURFACE FEATURE TYPES 


Integer Value Surface Feature 


No Elevation/No feature data 
Forest 
Urban 
Marsh 

Null(Elevation/No feature data) 
Water 

Heath(Waste land with shrubs) 
Open 


0 
1 
2 
3 
< 
5 
6 
7 





2. Node Characteristics 

In the network, nodes are used to represent actual terrain features on a map. 
Four characteristics are used to describe a node: an integer node identification number, 
three digit latitudinal coordinate, three digit longitudinal coordinate, and node type. 
The coordinates are simple map grid coordinates. The node type is entered as a single 
integer value. Eventually, each node has an altitude linearly interpolated from 100 
meter grid square data. Table 2 shows the integer values and the terrain feature it 
represents. 

The primary function of the node in the transportation network is to relate the 
junction of arcs in the data representation to physical locations on the map. In 
addition, they represent possible objectives in a combat mission. 

3. Arc Characteristics 

lhe arcs in the network represent feasible routes of maneuver from one 
location (node) to another. For a route to be considered feasible, it must provide for 
the maneuver of at least one column of tracked vehicles, wheeled vehicles, or 
dismounted troops. Three characteristics are used to represent an arc: head node 
identification, tail node identification, and arc type. These characteristics provide data 
for the algorithms to compute appropriate arc traversal times, distances and flow rates 


from node A to node B. 
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TABLE 2 
TYPE OF NODE 


Integer Value Terrain Feature 


City 

Village 

Road Junction 
ο. του 
Other 





The head node identification is the identification number of the node where 
the arc originates. The tail node identification is the identification number of the node 
where the arc terminates. The characteristics of the arc is an integer code for one of 


seven possible arc types as shown in Table 3. 







TABLE 3 
TYPE OF ARC 
Integer Value Terrain Feature 
1 Autobahn | 
2 Concrete or Asphalt Road 
3 Dirt Road 
4 Railroad 
5 Forest 
6 Open Country 
7 Bridge or Tunnel 
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4. Speed Data 
The speed for unit maneuver is a function of arc type, unit type and unit 
formation. The possible speed represents the speed on the specific arc for the specific 
unit for each unit formation. [t is assumed that the formation width of a single column 
is less than the width of the road itself. If the width of the unit formation is less than 
the width of the road itself, then arc speed is on-road speed. Otherwise, arc speed 15 
off-road speed, because on-road speed is greater than off-road speed with same arc 


type and unit type. 


C. STEPS IN PROCESS 
To translate 100 meter grid square data to attributes of arcs, the network 
structure is used. The network structure can exploit the mathematical nature of the . 
variables of the decision process by representing these variables as characteristics of the 
arcs and nodes in the network overlayed on the gridded terrain. [t is an objective in 
developing this network to adequately describe the terrain and combat effects on 
maneuverability necessary for the implementation. The process of translation and 
implementation from 100 meter grid square data to attributes of arcs is described for 
single arc geometry and later for multiple arc geometry and flow rates in the network. 
1. Single Arc Attributes 
a. Geometry 

Consider an arc which links node A to node B. To represent the 
coordinates of a node in three dimensions, the Cartesian coordinate (x, y, Z) system Is 
used. Each node has an altitude interpolated trom the 100 meter grid square data. 
From the coordinates of node A and node B, the total distance between node A and 
node B, and the angle of inclination of the arc is computed. The angle of inclination is 
determined by the slope between the x-axis and y-axis. The definition of angle of 
inclination of a line that crosses the x-axis is the smallest positive angle that the line 
makes with the positively directed x-axis [Ref. 3: p. 50]. 

If the angle of inclination is between 45 degrees and 135 degrees, then the 
X-axis is used as reference axis for calculation of distance between each point along the 
arc. Ifthe angle of inclination is greater than 135 degrees or less than 45 degrees, then 
the y-axis is used as the reference axis. Using this angle, coordinates of the cross 
points between a point on the arc and a point on the each 100 meter grid line for the 


reference axis is computed. After determining the coordinates of the points along the 
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arc, the distance between each pair of points is calculated step by step until the cross 
point reaches the tail node based on the distance formula for the plane. Figure 3.2 
shows an example of crossing points between the arc ( which is Node 1 to Node 2 ) 


and the x-axis 100 meter grid line along the arc. 
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Figure 3.2. Crossing Points. 
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b. Slope | | 

(1) Along The Arc. To determine the gradient pattern along the arc, the 
slope of the terrain between each pair of points along the arc is calculated. Coordinates 
and altitudes of each point along the arc and the difference of altitudes between each 
sequential pair of points are used to calculate the slopes between each pair of points 
along the arc. The coordinates of each point along the arc are computed by using the 
the point-slope equation of the line. As a measures of terrain steepness along the arc, 
three factors are computed: 


e Slope Change - sum of the absolute value of slopes between each pair of 
adjacent sub-arcs. 


e Total Slope Change - sum of slope change along the arc. 


e Average Slope Change - total slope change divided by the number of pair of 
sub-arcs along the arc. 


The slope of the surface along the arc is related to the feasible routes 
of maneuver for each type of unit from one node to another. Some arcs provide for 
the maneuver of dismounted troops, but tracked and wheeled vehicle movement тау 
be infeasible. Some arcs provide for the movement of tracked vehicles to climb steep 
slopes but wheeled vehicle movement is not possible. The maneuverability on a given 
slope depends on unit type ( tracked vehicles, wheeled vehicles or dismounted troops ) 
and soil strength which is measured by rating cone index points: the higher the index 
numbers, the greater the soil strength. The maximum negotiable slope can be used after 
obtaining soil strength data for arcs as shown in Appendix G. [Ref. 4: p. 5-19] 

Table 4 shows an example of general information and slopes along the 
arc from Node 1 to Node 2 as shown in Figure 3.2. XM represents the slope between 
starting node and terminal node from the coordinates of the x, y plane on the arc. 
ANGLE represents the angle of inclination of the arc. ACT.DIST is actual distance 
between two nodes along the surface, but DISTANCE is the distance on the plane or 
on the map. It is assumed that the distance of the arc is DISTANCE. AVERAGE 
SLOPE represents the difference between altitude of Node | and Node 2 divided by the 
distance of the arc. 

(2) Perpendicular to Arc. After determining the slope of the terrain along 
the arc, the slopes of the terrain along lines perpendicular to the arc are computed. The 
starting point is the cross point between the x-axis grid line and the arc, if the angle of 
inclination 1s greater than 45 degrees and less than 135 degrees. Otherwise, the starting 


point is the cross point between the y-axis grid line and the arc. 
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TABLE 4 
SLOPE AND DISTANCE 


NODE 2 NODE 1 ХМ ANGLE DISTANCE ACT. DIST 
2000 9000 8200 9800 -1.00 135.00 1131.37 1131.95 
POINTS ALONG ARC 
POSITION COORDINATE( X : Y) ALTITUDE 
POINT 1 3000.0 9000.0. 367.0 
POINT 2 8900.0 9100.0 267.0 
POINT 3 8800.0 9200.0 367.0 
POINT 4 8700.0 9300.0 372.0 
POINT 5 8600.0 9400.0 371.0 
POINT 6 8500.0 9500.0- 367.0 
POINT 7 8400.0 9600.0 361.0 
POINT 8 8300.0 9700.0 | 354.0 
POINT 9 ( 8200.0 9800.0 ) 348.0 
POSITION DELTA H DISTANCE SLOPE ACTUAL DIST 
POINT - 2 0.0 йаа 0.000 141.4 
POINT 2- 3 0.0 14144 0.000 141.4 
POINT 3- 4 580 141.4 0.035 141.5 
POINT 4 - 5 -1.0 141.4 -0.007 141.4 
POINT 5 - 6 -4.0 141.4 -0.028 141.5 
POINT 6 - 7 -6.0 141.4 -0.042 141.5 
POINT 7 - 2 -7.0 141.4 -0.049 141.6 
POINT 8- 9 -6.0 141.4 -0.942 141.5 
POSITION SLOPE CHANGE 
POINT 1- 2 0.00 
POINT 2- 3 0.04 
POINT 3- 4 0.04 
POINT 4- 5 0.02 
POINT 5 - 6 0.01 
POINT 6 - 7 0.01 | 
POINT 7 - 8 0.01 i 
TOTAL SLOPE CHANGE : 0.1697 
AVERAGE SLOPE CHANGE 0.0212 
NODE 2: ALTITUDE 1 NODE 1 ALTITUDE 2 
367.0 348. 
AVERAGE SLOPE -0.0168 





The altitude of the starting point is calculated by linear interpolation 
from the nearest left intersection between the x-axis and y-axis to the nearest right 
intersection between the x-axis and y-axis line. In the following discussion, it 1s 
assumed that the angle is 135 degrees. Table 5 and Figure 3.3 are used to illustrate the 


following discussion. 
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The sequence of calculations is to determine the slope of left off-road 
and right off-road along the perpendiculars. On the left side the first distance is from 
the starting point to the first cross point between the perpendicular line and y-axis grid 
line. The second distance is from the same starting point as the first one to the second 
cross point, etc. The altitude of the cross point between the perpendicular line to the 
arc and the y-axis grid line is computed by the same interpolation procedure as for the 
Starting point. Using the first distance and aititude of the starting point and the first 
cross point between the perpendicular line and y-axis grid line, the first slope is 
calculated. The process is continued until terminated by the stopping rule described 
later. After computing each slope on the perpendicular line for the initial point ( which 
is the head node ), slopes along the perpendicular for subsequent points along the arc 
are computed, terminating with the taii node. After determining the slope on the left 
side, the right side slopes are computed using the same procedures. Before calculating 
the slope on the perpendicular line to the arc, coordinates of each cross point between | 
the perpendicular hne to the arc and the y-axis grid line are determined. The 
coordinates of each point on the perpendicular line to the arc are computed using the 
same point-slope equation for the coordinates along the arc. Also, the difference in 
altitude of each pair of points is computed. 

lhe slopes on the perpendicular lines are related to the feasible arc 
width for movement from left boundary to right boundary. The arc width measured 
along the perpendicular line to the arc is considered for the trafficable path for each 
specific type of unit. Table 5 shows an example of slopes used to determine the 
trafficable width of an arc for the movement of vehicle units. For example, POINT 4 - 
| represents -the first (left) cross point between the y-axis 100 meter grid line and the 
line perpendicuiar to the arc at the fourth point from the starting node. POINT 4 - 5 
is the boundary point to left of the arc, POINT 4 - 6 1s the first point to the right, and 
POINT 4 - 11 is the boundary point to right of the arc as shown in Figure 3.3. 
SLOPE represents DELTA H divided by DISTANCE. 

(3) Arc Width. Two types of stopping rules are used for determining the 
arc width: a code change rule and a siope change rule. Briefly stated, an arc width 
boundary is established along a line perpendicular to an arc when the terrain code 
changes from type 1 to type j as described later. Also, a boundary may be established 


because the slope change along the perpendicular exceeds an input threshold value. 


24 


ТАВГЕ 5 
SLOPE ON PERPENDICULAR LINE TO ARC 


SOME POINTS OFF ARC (FROM EACH POINT ON THE ARC) 


COORDINATE ( X : Y) ALTITUDE DELTA H DISTANCE SLOPE 
PORNT 1- 1 8900 8900 374.00 7.0 141.42 0.05 
Joni Ι- 2 9100 9100 364.00 230 141.42 -0.02 
POTNIT 2- 1 8800 9000 374.00 7.0 141.42 0.05 
BOINI 2- 2 9000 9200 203.00 -164.0 141.42 О 
SOUNT 3- 1 8700 9100 374.00 7.0 141.42 0.05 
τ... 3- 2 8900 9300 366.00 =) 141.42 -0.01 
POINT 4- 1 8600 9200 371.00 το 141.42 -0.01 
BOINI 4- 2 3500 9100 384.00 12:0 282.84 0.04 
POINT 4- 3 8400 9000 404.00 32.0 424.26 0.08 
РОІМТ 4- 4 8300 8900 443.00 710 565.69 0 15 
POINT 4- 5 8200 8800 499.00 1270 707.11 0.18 
РОІМТ 4- 6 8800 9400 367.00 25-1 141.42 -0.04 
РОІМТ 4- 7 8900 9500 358.00 -14.0 282.84 ο 
POINT 4- 8 9000 9600 354.00 Δι 424.26 -0.04 
POINT 4- 9 9100 9700 35100 =21.0 Boon oe -0.04 
ӘСЕМ 4-10 9200 9800 348.00 -24.0 πο ЕЕ 204105 
РОІМТ 4-11 9300 9900 348.00 -24.0 848.53 -0.03 
ΕΠΙ 5-1 8500 9300 380.00 270 141.42 0.06 
POINT 5- 2 8400 9200 381.00 10.0 282.84 0.04 
БОЛИТ 5- 3 8300 9100 400.00 25170 424.26 0.07 
POINT 5- 4 8200 9000 440.00 69.0 565769 0.12 
BOL 5- 5 8100 8900 482.00 LIII ате O.16 
POINT 5- 6 8000 8800 499.00 128.0 848.53 0.15 
POINT 5- 7 7900 8700 499.00 128.0 83: 0:15 
ΕΙΝ 5= 8 7800 8600 492.00 121.0 зс DTE 
POTNI 5- 9 8700 9500 361.00 -10.0 141.42 ο 07 
НЕ 5-70 8800 9600 354.00 male ai) 282.84 -0.06 
POCNI 5-11 8900 9700 351.00 - 20.0 424.26 Zub 
(ουν 5-12 9000 9800 351.00 - 20 "0 295,099 -0.04 
БӘРІ 5-13 9100 9900 348.00 523.0 ο 11 -0.03 
РОТ 6- I 8400 9400 374.00 7,0 141.42 0.05 
POLINI бе 2 8300 9300 394.00 27.0 282.84 0.10 
POLNI o- 6S 8200 9200 398700 31.0 424.26 0.07 
POINT 6- 4 8100 9100 420.00 53.0 565.69 0.09 
БЕЛЕ” 6- 5 8000 9000 453.00 86.0 μα 212 
БЕП 6- 6 7900 8300 466.00 9910 848.53 0.12 
БОШКА $99 7800 8800 469.00 102.0 36595 0:10 
εως 6- 8 7700 8700 469.00 102.0 τοι η 0.09 
FOES 6- 9 8600 9600 3358700 -9m8 141.42 EU 
EOENT 6-10 8700 9700 351.00 EGO 282.84 -0.06 
ROINE 6-11 8800 9800 348.00 ZI О 424.26 -0.04 
ΠΠ, 6-12 8900 9900 348.00 το ο 565.69 -0.03 
rou 7- 1 8300 9500 367.00 6.0 141.42 0.04 
ROINE 7- 2 8200 9400 381.00 20.0 282.84 0.07 
τον 7- 3 8100 9300 384.00 Zone 424.26 0.05 
POINT 7- 4а 8000 9200 400.00 39.0 565,69 0.07 
РОРМТ =7- 5 7900 9100 404.00 43.0 707.11 0.06 
BOINI pi 6 7800 9000 407.00 46.0 848.53 0.05 
ΕΠΗ ΘΥ--ῃ 7700 8900 417.00 Sem SOS MD 0.06 
POINT 72 8 7600 8800 440.00 2-0 1231-37 0.07 
ΠΤΙ 7500 8700 449.00 88.0 1272.79 0.07 
POINT 7-10 8500 9700 τυ “10 Ὁ 141.42 2007 
Peni ет] 8600 9800 348.00 “15.0 282.84 -0 05 
ΙΕ 8700 9900 348.00 ze ῃ 424.26 70.03 
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TABLE 5 
SLOPE ON PERPENDICULAR LINE TO ARC (CONT D.) 


141.42 
282.84 
424.26 
369.09 
707.11 
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Тһе stopping rules for the arc width are a function of unit mission ( 
unit formation ) and unit type. When the untt formation 15 а single column, only the 
on-road portion of the arc is used. When the unit formation for movement is multiple 
column, both on-road and off-road attributes may be used. The arc width for a given 
perpendicular is the distance from the left side boundary to right side boundary, each 
established by one of the two rules. 

For example, in case of a vehicle unit, an open area ( code 7 ) and no 
feature data area ( code 4 ) may be considered feasible for routes of maneuver for off- 
road movement, with forest, urban, marsh, water, and heath areas considered 
infeasible. In the case of dismounted troops, however, forest, urban, heath and open 
area may be considered feasible for movement routes off-road, but water and marsh 
areas are infeasible. 

Consider the matrices shown in Tables 6 and 7. At the first cross point 
between the perpendicular line and y-axis grid line, it determines the code from the 
code of the nearest point on the y-axis grid line. If the value contained in the two 
dimensional array of the previous code and current code equals ‘no boundary, the 
process continues. If the value of previous code and current code is ‘boundary, it will 
stop. For example, a tracked vehicle can go from open area to open area but can not 
go from open area to water or marsh area. The dismounted unit can go from forest to 


forest, open or heath area, but will not go from forest to water or marsh area. That is, 
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Figure 3.3 Points on the Perpendicular. 


ING EC) equals 1, then the point is not a boundary, otherwise, the point is a 
boundary. Subscript i represents the previous code on the perpendicular line to the arc, 
and subscript j represents the current code on the same line as subscript 1. 

The slope, as well as terrain code change, could be a barrier to certain 
types of units. The slope of the perpendicular line to the arc affects the movement of 
units. The arc boundary due to slope occurs when the slope is less than a specified 


negative threshold or greater than a positive threshold. The threshold of slope is the 
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ТАВІЕ 6 
BOUNDARY FOR VEHICLE UNIT 
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TABLE 7 
BOUNDARY FOR DISMOUNTED TROOPS 
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maximum negotiabie slope for each unit type. The threshold for vehicle units may be 
less than that for dismounted troops. Eventually, arc soil strength data will be used to 


determine these thresholds as shown in Appendix G. 


28 


Each arc width on the perpendicular line is determined by 
consideration of the code change rule and slope change rule simultaneously. The arc 
width is the distance from the stopping point on the left side to the point on the right 
side of the arc. The line connecting each stopping point on the left ( right ) side is the 
arc’s left ( right ) boundary line as shown in Figure 3.4. 

2. Determine Arc Attributes 

Each arc width on the perpendicular line to the arc is determined by the code 
change and slope change rules previously discussed. The arc widths may be different 
for each perpendicular line along the arc. However, one unique arc width is required 
from a head node to a tail node, since it is difficult to change the unit formation during 
movement on an arc. Insertion of new nodes on the arc can be considered for long 
arcs when half of the arc is narrow but the other half of the arc is wide. 

The program determines the arc width by adding the left minimum distance, 
the right minimum distance, and road width itself. The minimum distance on the left 
side is determined Бу the shortest perpendicular distance from the arc to the left, and 
similarly for the right minimum distance. When the left and right minimum distances 
equal zero, the arc width is considered to be the road width itself since no off-road 
movement is possible for that unit type. The minimum distance determination of arc 
width is more realistic than consideration of some average of the distances. Тһе 
minimum distance has an direct effect on the flow rate of an arc. Table 8 shows an 
example of the minimum distance off arc, arc width and those coordinates for the plot 
in Figure 3.3. 

3. Flow Rate 

Using the arc width and possible speed for a specified unit type and mission, 
the rate of flow is computed. Possible speed is a function of arc type, unit type and 
unit mission type. The arc type represents several conditions of the maneuver paths 
between nodes and is an important element for the speed of maneuver. There are 
several categories of arc types as shown in Table 3 . The unit types currently 
considered are tracked vehicles, wheeled vehicles and dismounted troops as shown in 
Table 18, Appendix F. The types of unit formations ( missions ) currently considered 
are multiple column and single column as shown in Table 19, Appendix F. 

The following equation is used to determine the flow rate on an arc and 15 


measured in terms of battalions per hour for each unit type [Ref. 5: p. 20]. 


FR, 7 (SP, X W)/(DW,, x DD,,) (eqn 3.1) 
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Figure 3.4 Arc Width. 
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TABLE 8 
MINIMUM DISTANCE OFF ARC 


MINIMUM DISTANCE OFF ARC 


ШЕК iN DLSTANCE : 141.4 
RIGHT MIN DISTANCE š 141.4 
ROAD WIDTH : 670 
ARC WIDTH š 288.8 
POINT POINT 2 
SLEET Ү LEFT X RIGHT ІБН 
INITIAC POINT : 8900.0 8900.0 909590 9100.0 
ЕМІ) POINT =: 8100.0 3700.0 92950) 290070 


The equation is used for ail arcs, j, to compute the rate of flow, FRj, based on 
the speed of a battalion sized unit. SP, Eun width of the arc, Wi the doctrinal width 
of the battalion, DW,,_, and its doctrinal depth, DD,,. The doctrinal width and depth 
of the battalion depend on the unit type. Since the doctrinal width and depth of a 
battalion of dismounted troops is different from those of a tank battalion, different 
values were assigned for these parameters as a function of unit type. Table 9 shows an 
example of flow rate and maximum number of elements on the arc width for Arc 1 


(from Node 1 to Node 2) as shown in Figure 3.2. 


TABLE 9 
FLOW RATE 


ARC ТҮРЕ : 2( Concrete & Asphalt Road ) 
ОНЫТ ТУРЕ MISSION TYPE FLOW RATE МАХ + ELEMENT 
2.41 


Unit Туре! : Tracked Vehicle Unit 

Unit Type 2 : Wheeled Vehicle Unit | 
Mission Type 1 : Multiple Column Formation 
Mission Type 2 : Single Column Formation 
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= The maximum number of elements which can move in multiple column 
formation on the arc width is computed. The number is determined from the arc width 
and the distance between elements ( i.e., the maximum number of elements equal the 
arc width divide by the specified distance between elements for each formation type ). 
4, Arc Traversal Time 

For purposes of this planning model. it is not necessary to develop extremely 
accurate representations of attacking force movement. General estimates of expected 
movement rates based on unit tvpe and the characteristics of the terrain being crossed 
are sufficient. For simplicity, arc traversal time is the arc length divided by arc speed 
for the attacking force unit. 

Assume the formation width of a single column is less than the width of road 
itself. If the width of the unit formation is less than the width of the road itself, then 
arc speed is on-road speed. Otherwise, arc speed is off-road speed. On-road speed is 
assumed to be larger than off-road speed. Arc speeds, in terms of on-road and off-road 


speeds, are provided in Appendix A. 
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IV. NETWORK IMPLEMENTATION 


A. INTRODUCTION 

Methodology for determination of minimum time and minimum distance paths 
using a network representation of terrain is presented in this chapter. This chapter 15 
related to the computer program for the Cartesian space network in Appendix D. The 
Cartesian space network is used to find a shortest path between two given nodes on an 
undirected graph represented by a list of arcs. Input data for this network program 1s 


the output of the computer program of single arc attributes given in Appendix C. 
΄ 


В. DATA FILE 

The network representation of the maneuver area is placed on two disk files 
which are the output files from the computer program for single arc attributes: one for 
vehicle units and one for dismounted troops. Eight characteristics are used to represent 
an arc: arc identification, tail node identification, head node identification, traversal 
time, flow rate, length, width, and arc speed. These characteristics provide sufficient 
data for the algorithms to determine the minimum time or distance path from starting 
node to terminal node. The time for the arc is an integer representing the value of 
traversal time for the movement of each unit. The flow rate is an integer value 
representing flow rate for arc, |, times one hundred ( F К, x 100). The integer value 15 
used for computational reasons in the shortest path algorithm. The length is the 
distance between head node and tail node in kilometers. The width of the arc is the 
width of road itself plus the left and right off-road widths in kilometers for which the 
terrain will support movement. The arc speed represents the average speed for each 
unit along the arc. 

Two kinds of width determination rules are used for the input data files from the 
single arc attributes program. The width determination ruie is used in the following 
situation. Recall from Chapter III that either a terrain code or slope change may 
terminate the search for the arc width on each side of the arc. If the first point on a 
perpendicular to an arc terminates the search, the width determination rule is invoked. 
For example, suppose the distance to the first point is 100 meters. If rule 1 is used, the 
off-road width is assigned to be 100 meters for that perpendicular. If it 1s desired to 


consider off-road width to be zero, then Rule 2 1s used. 
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For example, if the unit was a convoy of wheeled vehicles, then Rule 2 would 
likely be selected. For tracked vehicles or dismounted troops, it may be appropriate to 
select Rule 1. Table 10 shows an example of the data for the network representation 


of vehicle unit from width determination rule 1. 


TABLE 10 
DATA FOR VEHICLE UNIT - RULE 1 


TAIL HEAD TIME FLOW Стар НЕЕ SPEED 
Unit* Unit* Km/Hr 
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Unit* : Hour x 100 
Unit** :( Battalions | Hour) X 100 
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TABLE 10 
DATA FOR VEHICLE UNIT - RULE 1 (CONT’D.) 
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C. DETERMINING THE PARAMETERS 

The program requires several parameters for determining the shortest distance or 
time path. Changeable parameters are unit type, unit formation, starting point and 
terminal point, arc travel time with obstacles, and whether minimum time or distance 1s 
to be calculated. During the program runs, the parameters are provided by user 


interaction with the program. 
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1. Unit and Path 

The path can be selected in two ways: minimum time or minimum distance. 
Minimum time determines the least time route through a network from a user selected 
starting node, ISTART, to a terminal node, LAST. Minimum distance determines the 
shortest distance path between two selected nodes. 

The unit type is vehicle unit or dismounted troops based on a battalion sized 
unit. Vehicle unit 1s assumed to be tracked vehicles for the examples described later. 
When the unit is determined, the program reads the data file appropriate for the 
selected unit 

2. Unit Formation 

The unit formation is described as either a single or multiple column 
formation. For dismounted troops, a single column formation is defined as two 
columns, with muitiple column formation being more than two columns. 

The width of the unit formation is the number of columns times the input 
spacing between columns. The depth of the unit formation is the number of rows times 
the input spacing between rows. If the width of the formation is larger than the arc 
width, then the program assigns a “large” number to the travel time, which is the arc | 
cost in the shortest path aigorithm, Otherwise, the program uses the original value for 
time previously determined as the arc cost. | 

3. Obstacles | 

Obstacies ( e.g., minefields ) can change the travel time and the shortest path 
in the program. The program requires information about obstacles on the arcs. If 
there 1s an obstacle on an arc, then the travel time is appropriately increased. For the 
purpose of examples in this. thesis, the original travel time is multiplied by four when 


obstacles are present. 


D. SHORTEST PATH ALGORITHM 

The determination of the minimum time and distance path through the sub- 
network is one of a class of many such problems often referred to as “shortest path” 
network problem. The general shortest path problem is to determine the least cost 
route through a network starting at node ISTART and ending at node LAST. The 
cost of a route 1s some function of the characteristics of the arcs and nodes that make 
up the route from node ISTART to node LAST. In the case of this minimum time 


path problem, the cost of traversing an arc is the amount of time it takes a unit to 
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traverse the arc. The cost of the route is the sum of the costs of all the arcs in the path 
from the supply node to the terminal node. 

Before selecting a method of finding the shortest path, the network must be 
examined to insure that there is, in fact, a solution in all cases. Because no arc in the 
network will have a negative cost associated with it, the addition of another arc to any 
path will never reduce the total time traveled. Therefore, negative cycles cannot exist in 
the network. The only arcs which have a zero cost are those leaving the supply node 
and entering the terminal node. Furthermore, boundaries are selected to run parallel to 
avenues of approach through a sector, so it is valid to assume that there will be at least 
one set of connected nodes which extends laterally through the sector. Thus, the 
method of construction of the sub-networks insures there will always exist at least one 
path between supply and terminal nodes. Therefore, there must be a solution to the 
shortest path problem, and there will not be any cycles in a minimum path( no node 
wil be visited twice). [If the assumption that boundaries are drawn which include a 
connected path from the start to the terminal node 15 invalid, a solution may not exist 
to the problem. 

Because of the many applications for this class of problems, there are many 
algorithms available for its solution. Most of these algorithms consist of two 
procedures: a label correcting procedure and a search procedure. 

In the label correcting procedure each node is initially assigned an infinite cost 
with the exception of the starting node, s which is given a cost of zero. Letting c(u) be 
the cost assigned to node u, d(u,v) be the cost to traverse the arc from node u to node 
v, and pred(u) be the node previous to u in a path, then the following rule is applied to 


change the costs associated with node v: 


If c(u) + d(u,v) < c(v) (eqn 4.1) 
then c(v) = c(u) + d(u,v) and pred(v) = u. 


Though inefficient if applied indiscriminately, if this rule is continuously applied 
until no cases of Equation 4.1 being true can be found, the chain of pred(v) to pred(u) 
= s will describe the minimum path from every node in the network to the starting 
modes. |Кег 6: p. 61] 

The search procedure is used to decide in which order the nodes are to be 


scanned to apply the labeling rule. This step is where most algorithms differ and also 
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where the efficiency of the algorithm is achieved. In most cases one of three search 
techniques are used: Dijkstra's algorithm, depth first search, and breadth first search, 
but many hybrid techniques are also available. Dijkstra’s algorithm gives priority of 
search to the adjacent node which is the shortest distance away. Depth first search 
gives priority to the most recent node searched. Breadth first search gives priority to 
the oldest node searched. The Dijkstra's algorithm was selected for use in the program 
because all arc lengths are at least greater than zero and are not the same. 

As the name implies, the shortest path algorithm solves arc cost in terms of the 
minimum travel time and distance for an attacker across an arc. For simplicity, the 
attacker force is treated as a rectangle for all computations and travel time is the arc 
length divided by average arc speed of the attacking force. The former is an arc 
attribute, while the latter is determined from unit movement capabilities. The algorithm 
is described below. 

Input: Cartesian Space network G(V,E), nodes V, arcs E with their 
characteristics as described in Section B. 
output: Shortest path distances or times( d(v) ) from starting node, s, 
to terminal node, t. 
Step 1. Initialization: | 
a) set all labels to a very large value. 
d(v) = infinity 


b) set cost of starting node to zero. 


d(s) = 0 
c) place starting node in the Set S. 
Set S = {s} 


b) remove the Set S from Set V, the set of nodes of G. 
Set T = Set V - {s} 
Step 2. For each node v adjacent to the node u (s.t. u e S, ve T). 
(If ао) + u,v} < (У) 
‘then (У) = d(u) + (u,v) } 
Step 3. Determination of Set S and Set T. 
a) determine V hin 
V nin > argmin (d(v); 
b) save the value of cost (d(v)) in the Set S. 
бенокил: 
c) remove the Set S from the set of nodes of G ( Set T ). 


38 


Set T = Set T - (5) 

d) If all arcs in the graph have not been analyzed, go to Step 2. 
Step 5. Return Set S to the minimum time path solution routine. Stop. 

(u,v) is the length of the arc between node u and adjacent node v. 

Set S is the set of nodes such that d(v) gives the true shortest 

path length from starting node s to some node v. 

Set T is Set V - Set S 

Set V is the set of nodes of undirected graph G. [Ref. 7: p. 129] 

When Set T is empty, each node in the network will have a cost assigned which 
represents the minimum cost for all possible paths from the node to the starting node.s. 
The minimum path from any node to the start node can then be found by tracing the 
chain of predecessor values back to the start node. When the algorithm is completed, 
the minimum time or distance path is defined by the chain of predecessor values from 
the supply node,s, to the terminal node. After determining the node number aiong the 
minimum path, the arc number is found for computing total travei time or total 
distance from the start node to the terminal node. The total travel time for the unit is 


computed as follows. 
Ет ә) + (UL / SB) (eqn 4.2) 


The equation is used for total travel time for a unit, T nip based on the total 


.) and unit length, UL, divided Бу 


агс 1 
average speed along the path, SP. The second term of Equation 4.2 represents 


travel time for an element( summation of T 


additional travel time for the entire unit to complete the path. 


The following equation is used to determine average speed along the path. 
a= > (dad x V )/ D | (eqn 4.3) 


The average speed SP is computed based on the summation of length of arc., d., 
multiplied by the speed on arc., V., and that quantity divided by the total path length, 
D. This algorithm finds the minimum time or distance path through the sector by 
using the unit formation. Any type of arc can be on this path and speeds used to 
determine travel time are adjusted for the type of arcs along the minimum time path. 
Table 11 shows an example of results for the vehicle unit from the network program 
given in Appendix D. 
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TABLE I 
OUTPUT FOR Εμ ΕΝ Rug cee 


ap аз аз ap ез ap == => => => => аз ap ез ез ер ер => ер ез ер «πὸ αν αν αν αν αν αν αν αν αν αν αν αν αν => ер => => => => ер ер => = = 


VEHICLE UNIT 
FORMATION : 50 ROWS 1 COLUMNS 
ROW SPACE 50 METERS 


NO MINEFIELD ARC 


SUM OF TRAVERSAL TIME : О HOUR 36 MINUTE 
TOTAL TRAVERSAL TIME : 0 HOUR 42 MINUTE 
NODE NUMBER ALONG MIN TIME PATH 

N ΝΌΡΕ 

1 1 

2 2 

3 3 

4, 4, 

5 5 

6 12 

7 16 

8 23 

9 29 

ШЕ 21 

ARC NUMBER ALONG MIN TIME PATH 

М АВС 

1 1 

2 3 

3 5 

4 4 

5 13 

60031 

τ; 

з 56 

9 68 
TOTAL РАТН LENGTH 14.60 KM 
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V. ANALYSIS 


A. INTRODUCTION 

This chapter describes the results of varying parameter values of the single arc 
attributes program and Cartesian Space Network program used in determining 
minimum distance or minimum time path. The results in this chapter are final outputs 
from the single arc attributes program through the Cartesian Space Network program. 
A 10 X 10 Km sector of European terrain is used for the analysis to demonstrate the 
algorithms. There are 34 nodes and 72 arcs in the transportation network. The start 
node is Node | and the terminal node is Node 34. The unit size is battalion level 
consisting of 50 tanks in a tracked vehicle unit and 1000 men in a dismounted troop 
unit. Minimum distance path is first discussed, followed by a comparative analysis of 


the minimum time path generated under various conditions. 


B. MINIMUM DISTANCE PATH 

This algorithm determines the path that provides the shortest distance through 
the network. The arc cost is defined by the distance between head node and tail node. 
The tactical significance of this algorithm would be the determination of the path 
through the network that provides the best route to reach the given point in terms of 
distance. A path of this nature might be used by raiding or reconnaissance elements of 
a dismounted troop unit. 

The minimum distance path is Path D through the network shown in Figure 5.1 
The distance of the path is 10.9 Km from the start node ( Node 1 ) to the terminal 
node ( Node 34 ). Figure 5.1 shows the minimum distance path and several minimum 


time paths described in the next section. 


C. MINIMUM TIME PATH 
The purpose of this section is to demonstrate the full capabilities. of the 
algorithms developed in terms of minimum time paths. The following parameters are 
considered in the analysis: 
ο Unit type - vehicles and dismounted troops. 
e Formation size - specified by number of rows and columns in the formation. 


e Width determination rule - specifies the method used to determine total arc 
width in the special case described in Chapter IV. 


e Obstacle ( minefield ) representation. 
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MIN DISTANCE AND TIME PATH 
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Figure 5.1 Min Distance and Time Path. 


42 


e Bridge representation. 
The results of several model runs for various parameter values are given in Tables 
12 and 13. The total travel time for the lead element and the unit, the path length, and 
the minimum time path from Figure 5.1 are presented for each case. The impact of the 
width determination rule selected is shown in Table 12 for vehicie units and in Table 13 


for dismounted troop units. 


TABLE 12 
MIN TIME PATH - VEHICLE UNIT 


Mine |Bridge| Travel Time! Length | Path 
Km 


y/n y/n jelement/unit 













А. А 
А. À 
Ы” А 
pe 
B 
B 
B 
B 
B 
3 
JO T2 14.6 À 
no feasible; route 
no feasible! route 
no feasible; route 
Ελ sox | 2 ves no 568} 13. | 14.2 | 8 
55... 2 2 yes no no ieasible| route 
ШӘ | 109 5 2 yes no no feasible| route 
E.4 > A 10 2 yes no no feasible| route 
36 43 14.2 B 









ves no feasible 
yes no feasible 
yes no feasible 
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1. Vehicle Unit 
The formation of vehicle units is either a single column or multiple column 
formation. The width of the formation for a single column 1s five meters, while that for 
the multiple column formation 1s 25 meters for each column. For example, the width of 
a 50 Х 1 formation is five meters, of a 25 X 2 formation is 25 meters, and of a lO X 5 
formation is 100 meters. Note from Table 12 that for vehicle units, Path A and Path B 
are selected for the minimum time path under various situations, because these paths 
do not have arcs which are in a forest or marsh area. 
a. No Minefieid 
For the no minefield and no bridge situation, the minimum time: path 1s 
Path A in the two cases which use width determination rules 1 and 2 ( see Cases А.1, 
A.2, A.3, and D.1 ). The path from the results of Rule l yields an optimal route for 
movement until the formation size is 10 X 5 ( ie., formation width is 100 meters in 
which case no feasible route exists( see Case A.4 in Table 12 }). 
.. For Rule 2, a path exists only for the single column formation ( Case D.1 ). 
For wider formations, there 1s no feasible route for movement ( see Cases D.2, D.3, 
D.4 in Table 12 ). The distance of the minimum time path is 14.6 Km and travel time 
for an element is 36 minutes, as shown in Table 12 . Note that the unit travel times 
vary as a function of formation size. 
b. Minefield 
The purpose of runs labeled Cases B and E in Table 12 is to demonstrate 
the effect of insertion of obstacles on the network. For these runs, a minefield is 
positioned on Arc 8 which is between node 5 and node 12 in Figure 5.1. Path B 
bypasses the arc which has a minefield for both Rule | and 2. The reason 1s that, for 
these runs, arc cost is four times larger when a minefield is on the arc as compared to 
the no minefield case ( 1.e., the travel time is four times larger for a minefield arc ). 
Other than a change from Path A to Path B, the effects of the width 
determination rule and formation size are similar to the no minefield case. 
c. Bridges 
Bridges tend to constrain movement of units to only the on-road portion of 
the arc. For these runs ( shown as Cases C and F in Table 12 ), a bridge was inserted 
on arc 8 instead of a minefield. Even though Rule 1 ( Case C runs ) are included, only 
Rule 2 would be applied in the real situation. The minimum time path bypasses the 
arc which has the bridge for both rules. The width of the bridge is assumed to be four 
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meters. Because the width of the bridge is less than the width of the formation in all 
cases, the minimum time path algorithm does not include the bridge arc on the optimal 
path. The results for this case are the same as for the minefield case. 
2. Dismounted Troop Units | 

The formations for dismounted troop units are considered as either one double 
column ( six meters wide ) or multiple column ( five meters for each column ). For 
example, a 500 X 2 formation is six meters wide, while a 50 X 20 formation is 100 
meters wide. For a comparison of vehicle versus troop units with the same formation 
widths, consider Cases A.2, G.2 and D.2, F.2 in Tables 12 and 13. Note that under 
Rule 1, Path C is optimal for troop units ( Case G.2 ) while Path B 15 optimal for 
vehicle units ( Case A.2 ). For the wider formations ( Cases G.4 through G.6 ) the 
optimal path shifts to Path B under Rule 1. For the tighter arc width restrictions of 
Rule 2, Path B is optimal in all cases for dismounted troop units. 

Finally, the effect of long units for travel time is evident from Case G.1. Note 
that an additional 44 minutes is required for the entire unit to complete Path C. 

This chapter has demonstrated various aspects of the algorithms developed in 
this thesis. Except for those areas described in Chapter VI, the model is ready for use 
in a full production mode for any terrain aréas which have digitized terrain data and 


transportation networks in the form described in Chapter ITI. 
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VI. SUMMARY AND FUTURE RESEARCH AREAS 


It is essenual that any combat simulations of the Airland battle provide a 
realistic representation of terrain and environment. Previous approaches have typically 
utilized either hex or grid squares within the model to describe terrain conditions. In 
some cases network overlays were used for convoy movements, but no models exist 
which capture essential terrain characteristics without the internal use of some grid 
pattern. 

In order to take advantage of the various algorithms available using networks 
without the need for gridded terrain representation, it is necessary to describe essential 
terrain parameters as arc and node attributes. A family of algorithms were developed 
and demonstrated in this thesis which accomplishes this objective. These algorithms 
can be implemented for any terrain area for which appropriate data exists. The 
resulting network is then usable in a model such as ALARM for both the planning and 
execution phases of combat representation. 

Areas of future research include the development of a maximum flow path 
algorithm. This area requires a significant effort in determining appropriate arc costs 
and development of maximum flow methodology for an undirected network. [а 
addition, several refinements are possible to represent other aspects of terrain 
characteristics not explicitly contained in the current model. Finally, additional test 
runs of the model with different terrain areas, unit types, combat formations,and 


boundary ruies need to be made. 
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APPENDIX A 
ARC WIDTH AND SPEED 


The detailed modei output for the network shown in Figure 5.1 is given in this 


appendix. 


TABLE 14 
VEHICLE UNIT - RULE 1 


Arc Node Road Off-Road | Off-Road] Total Speed 
H m Width Left Right Arc Width| On off 


1 ] 2 6 141 141 288 5 

2 1 8 6 100 7 113 40 25 

3 ο. 5 109 12 127 40 25 

4 2 9 2 4 4 10 5 5 

5 2 10 2 108 10 120 5 5 

5 з 4 5 111 0 117 40 25 

7 3 10 2 | 1 120 123 5 5 

8 4 5 5 142 2-56 174 40 25 

9 4 10 5 et ? 123 40 25 
10 ШЕШІ 5 6 104 115 300 20 
11 5 6 5 133 . 83 222 40 25 
12 5 11 4 119 11 222134 308 20 
13 5 12 6 9 108 123 40 25 
14 6 12 5 101 11 118 40 25 
15 7 6 4, 129 22 155 30 2 
16 ТР 4 123 0 127 30 20 
17 ЕСЕ 4 101 0 106 30% 20 
18 8 9 5 102 0 107 30 20 
19 8 18 2 5 100 108 5 5 
20 9 10 2 17» 56 170 5 Ε 
21 9 14 2 101 13 115 5 5 
22 9 18 2 137 9 148 5 5 
23 9 19 2 5 108 115 5 5. 
24 HORE 6 ШІ 100 100 206 40 25 
25 10 14 5 115 uo 30 20 
26 10 15 2 8 113 123 5 5 
27 ШІ 12 5 100 105 30 20 
28 11 15 5 105 110 40 25 
29 11 16 26 15 119 138 40 25 
30 128 13 4 141 0 145 3 20 


Unit of on(off)-road speed : Km | Hr 


Units of several widths : Meters 
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TABLE 14 
VEHICLE UNIT - RULE 1 (CONT’D.) 


6 101 7 114 
2 116 7 125 
2 114 124 
6 100 38 144 
2 122 16 140 
5 100 0 105 
+ 100 10 113 
2 103 0 105 
6 100 100 206 
2 128 0 130 
6 114 10 130 
2 108 0 110 
2 138 7 147 
2 100 107 
6 141 141 289 
2 228 62 292 
6 130 8 143 
6. 109 49 164 
6 107 15 126 
2 16 128 147 
2 12 109 124 
5 103 26 134 
5 14 134 153 
6 117 23 146 
6 115 13 СЕ 
6 109 ЗТ 146 
5 I3 101 118 
4 56 112 172 
2 Э 117 129 
2 119 10 131 
2 100 100 202 
5 1112 56 173 
6 6 101 113 
2 153 7 147 
2 141 141 255 
2 100 100 202 
2 108 0 110 
6 | 0 104 110 
6 

6 

+ 

2 
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TABLE 15 
VEHICLE UNIT - RULE 2 


Arc Node Road {Off-Road | Off-Road, Total Speed 
H T Width Left Right Arc Width! On ОБЕ 
40 25 


1 2 6 0 0 6 
-2 1 8 Š 0 0 6 40 25 
3 Ж 3 5 0 0 6 40 25 
4 2 9 2 0 0 2 5 5 
5 2 10 2 0 0 2 SÉ 5 
6 c M 5 0 0 6 40 25 
7 3 10 2 0 0 2 5 5 
8 4 5 6 0 0 6 40 25 
9 4 10 6 0 0 5 40 25 
10 4 11 5 0 0 5 30 20 
11 5 6 6 0 0 6 40 25 
12 Š 11 4 0 0 4 30 20 
13 5 12 6 0 0 6 40 25 
14 6 12 6 0 0 5 40 25 
15 7M 6 4 0 0 4 30 20 
16 7 12 4 0 0 4 30- 20 
17 ἽΝ 15 4 0 0 4 30 20 
18 8 9 5 0 0 5 30 20} 
19 8 18 2 0 0 2 5 5 
20 9 10 2 0 0 2 5 5 
21 9 14 2 0 0 2 5 5 
22 9 18 2 0 0 2 5 5 
23 9 19 2 0 0 2 5 5 | 
24712 ій ІТ 6 0 0 6 40 25 
25 106 та 5 0 0 5 30 20 
26 10 15 2 0 0 2 5 
27 111 12 5 0 0 5 30 26 
28" ||| 11. 15 6 0 0 6 40 25. 
29 11 16 6 0 0 6 40 25 
30 П 15 4 0 0 4 30 20 
31 12 16 5 0 0 6 40 25 
32 13 16 2 | 0 0 2 5 5 
33 ise ig 2 | 0 0 2 5 5 
34 14 15 6 | 0 0 6 40 25 
[ 35 14 19 | 2 | 0 0 2 oe τς 
36 158 16: || 5 0 0 5 30 20 
37 158 21 4 0 0 4 30 20 
38 16 17 2 0 0 2 5 5 
39 16 23 6 0 0 6 40 25 
40 16 24 2 0 0 2 5 
| 41 ІІ απ 6 0 0 6 40 25 
42 18 19 2 0 0 2 5 5 
43 18 26 2 0 0 2 5 5 
АН | 15. αι 2 0 0 2 5 5 
45 19 20 5 0 0 6 40 25 
46 19 26 2 0 0 2 5 
47 20 21 6 0 0 6 40 25 
48 | 20 22 6 0 0 6 40 25 
49 | 20 27 6 0 0 5 40 25 
50 | 20 28 2 0 0 2 5 5 
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ТАВГЕ 15 
VEHICLE UNIT - RULE 2 (CONT'D.) 


2 
5 
5 
6 
6 
6 
ч 
4 
2 
2 
2 
5 
6 
2 
2 
2 
2 
6 
6 
6 
4 
2 
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TABLE 16 
DISMOUNTED TROOP UNIT - RULE I (CONT D.) 


2 
5 
5 
6 
6 
6 
5 
4 
2 
2 
2 
5 
6 
2 
2 
2 
2 
6 
6 
6 
4 
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TABLE 17 


DISMOUNTED TROOP UNIT - RULE 2 
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ТАВГЕ 17 
DISMOUNTED TROOP UNIT - RULE 2 (CONT D.) 
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APPENDIX B 
DOCUMENTATION FOR COMPUTER PROGRAM 


SINGLE ARC ATTRIBUTES 
ЖЖЖ Ж ЖЕЖ DOCUMENTATION ote oie ots oie ois fe ois oie aie ois ЖЕ 
This FORTRAN program determines the attributes for individual arc from the 


100 meter grid square data, node characteristics and arc characteristics of the network 
overlayed on the gridded terrain. 
Assumption : 


ә Ifthe angle of inclination is between 45 degrees and 135 degrees, then use X- 
axis. m 


ο If the angle of inclination is greater than 135 degrees or less than 45 degrees, 
then use the Y-axis. 


FUNCTION XCORD: Compute X grid coordinates. 
FUNCTION YCORD : Compute Y grid coordinates. 
SUBROUTINE INIT : Set the initial conditions. 
SUBROUTINE INPUTI: Read the data from data 1, data 2, and data 3. 
SUBROUTINE INPUT2: Read the data about the arc. 
SUBROUTINE DST : Compute distance along arc. 
SUBROUTINE SLOP! : Compute slopes between each pair of point along arc. 
SUBROUTINE SLOP2 : Compute slopes between each pair of point off arc. 
SUBROUTINE MIND : Compute minimum boundary distance from the arc. 
SUBROUTINE FLOW : Compute flow rate for each arc. 
SUBROUTINE TIME : Compute min traversal time for each arc. 
SUBROUTINE PRT : Print the results of calculation. 
Ж жок жок эй >к жй Жок окси VARIABLE DEFINITION Жж эк э эк Р дәй жй жок дж 
CASE : option of military unit type 

CASE 1; vehicle unit 

CASE 2; dismounted troops 
PRINT : option of print 

PRINT 1; print all results of calculation. 

PRINT 2; print the minimum distance off arc. 

PRINT 3; print the input data for checking. 

PRINT 4; print the coordinates of boundary line. 

PRINT 5; print the data for network. 
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UNITF : type of formation 
UNITF 1; multiple column formation 
UNITF 2; a single column formation 
THRESHOLD of slope for boundary line 
THRESI : threshold for positive slope ` 
THRES2 : threshold for negative slope 
ARCTP : arc type from a head node to a tail node. 
ACTD  : actual surface distance of each pair of points. 
AD  :total actual surface distance along the arc. 
ADIST : distance between head node and tail node (in kilometers). 
ALTD  .: altitude of each point ( raw data ). 
ALTI .: altitude of the head node on the arc. 
ALT2 : altitude of the tail node on the arc. 
ANGLE : angle between arc and X-axis (in degree). 
AREAL : real variable by using array. 
ASLPCH : average slope change along the arc. 
ATIME : minimum traversal time for each arc. 
AVGSLP : average slope along the arc. 
DELTAX : distance between node 1 and node 2 along X-axis. 
DELTAY : distance between node 1 and node 2 along Y-axis. 
DIST : distance between node 1 and node 2 on the map. 
DISTA : distance between each pair of point along arc. 
DISTB : distance between each pair of point off arc. 
DHT _ : difference between height of each pair of point along arc. 
DHT2 .: difference between height of each pair of point off arc. 
DSTN .: distance of each element in line (in meters). 
DW _ : doctrinal width of battalion during deployment. 
DD _ : doctrinal depth of battalion during deployment. 
FORWV : minimum formation width of vehicle unit. 
FORWD : minimum formation width of dismounted troops. 
FR _ : flow rate of the arc (battalion / hour). 
HT _ : height of intercept between arc and Y-axis along arc. 
HT2 : height of intercept between arc and X-axis off arc. 
INTS : single integer variable. 


INTA  : integer variable by using array. 
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IX,IY : crossing points between an arc and 100 meter X, Y grid line. 
KNODE : number of crossing point along the arc. 

KR _ : coefficient of route availability. 

LMIN : minimum distance of the left side off arc. 

MAXNO : max number of elements that move in line along the arc. 
МАКС : total number of arc in the data file. 

NCODE : characteristics of each point ( raw data ). 

NCODEX : dummy variable for characteristics of each point. 
NODE .: total number of node in the data file. 

NODEA : identification number of head node on each arc. 
NODEB : identification number of tail node on each arc. 

NSPCL : variable of using width determination rule 2 for boundary. 
NX  : dummy variable for altitude of each point ( raw data ). 
RMIN : minimum distance of the right side off arc. 

ROADW : width of road itself. 

SLOPEI : slope between each pair of point along arc. 

SLOPE?2 : slope between each pair of point off arc. 

SLPCH : slope change between each pair of point along the arc. 
SREAL : single real variable. 

SP  : movement speed for an element of unit. 

TSLPCH : total slope change along the arc. 

UTMX : U.T.M. grid coordinates of the point ( East - West ). 
UTMY : U.T.M. grid coordinates of the point ( South-North ). 
WIDTH : width of arc ( road width + off-road width ). 

XA,YA : intercept point between arc and Y-axis. 

XB,YB : intercept point between reference line and X-axis. 

XL,YL : X,Y coordinates of minimum distance of left side off arc. 
AR,YR : X,Y coordinates of minimum distance of right side off arc. 
ANODE : X coordinates of each node. 

YNODE : Y coordinates of each node. 

XM ~ : slope between X-axis and an arc. 


ҮМ  : slope between Y-axis and an arc. 
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CARTESIAN SPACE NETWORK 
ck ok dex Ἂς Ἂς * * κ x DOCUMENTATION ЖЖЖ oe ook ox k kkk kkk kk k 
This FORTRAN program determines the minimum time path and minimum 

distance path from the undirected graph by using a shortest path algorithm. Input 
data of this program is output data from the single arc attributes program. 

SUBROUTINE OPTION: Determine the initial options. 

SUBROUTINE INIT : Set the initial conditions. 

SUBROUTINE FORMTN: Specify the formation width and depth. 

SUBROUTINE SHORTP: Determine the minimum time or distance path. 

SUBROUTINE PRINT : Print the results of calculation. 


жжжжжжжжжжжж ж VARIABLE DEFINITION ж Жк жок ж жож ok o ok oe ok ok 


CASE : Option of Military Unit Type 
CASE 1; vehicle unit 
CASE 2; dismounted troop unit 
UNITE: Type of Formation 
UNITF 1; multiple column formation 
UNITF 2; a single column formation 
PARAM : Parameter for Shortest Path 
PARAM 1; minimum time path 
PARAM 2; minimum distance path 
BIG : an arbitrary large real number. 
DIST : distance between head node and tail node on the map. 
FLOW : flow rate of each arc. 
HEAD  : head node of each arc. 
INODE : tail node of each arc. 
[NTS : single integer variable. 
INTA :integer variable by using array. 
ISPATH : integer vector of length iN. 
ISTART : start node of the minimum time or distance path. 
IWORK?2 : integer vector of length N. 
IWORK23 : integer vector of length M. 
JNODE : head node of each arc. 
LARC .:arc number along the minimum parameter path. 
LAST: terminal node of the minimum time or distance path. 
LINK :arc number in the graph ( link(head node, tail node) ). 
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М : number of arcs ( edges ) in the graph. 

MINE_: are number which 15 on the minefield. 

N : number of nodes in the graph. 

NARC : total number of arc in the data file. 

NODE .:total number of node in the data file. 

NO  :arc number in the graph ( one dimensional array ). 

NP  :the integer NP has the value zero if a shortest path is found. 
between ISTART and LAST, otherwise it has the value one.. 

NUMP .:the number of nodes in the shortest path found. between. 
ISTART and LAST. 

PATHL : total-distance between ISTART and LAST. 

REALS : single real variable. 

REALA : real variable by using array. 

RESM : character variable RESM has the value Y if there is a mine’ 
field in the graph, otherwise it has the value N. 

TAIL : taillnode Oi cach aro: 

TIME : minimum traversal time of each arc. 

TPATH : total traversal time for the entire unit between ISTART 
and LAST. 

WIDTH : width of arc ( road width + off-road width ). 

WK4  :real vector of length N. 

XLEN :total value of parameter between ISTART and LAST. 
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APPENDIX C 
COMPUTER PROGRAMI FOR SINGLE ARC ATTRIBUTES 


This appendix contains the computer program used to support the algorithms in 
Chapter III. 
PROGRAM TSIM 
JULY 5 1987 (20:00) 
kkXAdxkkkk** VARIABLE DECLARATION ХХХХЯХААААЯЯ 


INTEGER  ALTD,ARCTP,CASE,CODE,PRINT, UTMX,UTMY,UNITF 
REAL LMIN 

CHARACTER  *1 BLANK 

DIMENSION ALTD(100,100) ,NCODE(100,100 NX (10) ,NCODEX(10) 
DIMENSION ME „УМОВЕ (90) , СОВЕ (0:7,0:7) 

DIMENSION DISTA 90) , YA(90 

DIMENSION HT (90) DHT (30), 156851 (60 

DIMENSION DISTB(100,100 р Ко, , YB(100,100) 
DIMENSION  HT2(100,100),DHT2(100,100),SLOPE2(100,100) 
DIMENSION МО 20), , SLPCH(90 
DIMENSION XL(10),YL(10), πα Жадо), 


* 
* 


1 
DIMENSION ΡΕ 5,2), ο), 
DIMENSION ACTD(90),ATIME(3 
DIMENSION ROADW(7 
COMMON / INTS / NSTEP,IX1,IY1,1IX2,1¥2,JA,JB,K,NARC,NCOUNT, 
P NODEA,NODÉB,Ji,J2,J3,J4,PRINT,UTMX,UTMY, 


= CASE, ARCTP, FR, AD, UNITF, NSPCL 
COMMON / INTA / NCODÉ, NO, ЖЕП. CODE, МАХМО 
COMMON / SREAL / XM, ANGLE , DIST, TSLPCH, ASLPCH,ALT1,ALTZ,AVGSLP, 

Š LMIN, RMIN, ЕГ YL, 3R, YR, ЕКТ, WIDTH, THRES1, THRES 2 
COMMON / AREAL / XA, YA, Н DHT, DISTA ,DISTB, SLPCH, SLOPE1, SLÓPE2, 

Ж БУВ ,HT2, DHT2, ЭР” απο. АТІМЕ, ROADW 


ШЕ СЕТ ТІНЕ INITIAL OPTIONS. 


WRITE(6,601) 

601 FORMAT ( ІІІ 

16 PRINT* ae | WHICH ТҮРЕ OF UNIT DO YOU DESIRE? ( 1 OR 2 )' 
PRINTA| 1 Па кв 
PRINT*, 1 * ЖІ 
PRINT*, | * 1. VEHICLE UNIT жі 
PRINT*; 1 * X! 
PRINT*, | ж 2. DISMOUNTED TROOPS ж! 
PRINI*, 1 ж жі 
PRINT ! Zç Zç Zç Zç Z Z Z Z Zç Zç Z < Zc Z Z Z Zç Je ke ke k Ze Z Z< Ze ZX Ze Zç Ze Ze Zç 1 
ο ж) 


А<Е i 
PRENT κα NOTE : FOR YOUR REFERENCE, CURRENT ANSWER IS', 
ASE 
IF ( "CASE κο. ο. CASES EQ: 2) THEN 


COLTO 17 
ELSE IF ( CASE .NE. 1 .OR. CASE .NE. 2 ) THEN 
PRINT* 
PRINT*, AAR ERROR : ENTER THE NUMBER 1 OR 2 ***' 
PRINT*, ' 
59 ТО 16 
END. IE 
17 WRITE(6,602) 
602 FORMAT M ὄν. 
PRINT*, |! eS TYPE OF FORMATION DO YOU DESIRE?' 
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18 
603 


604 


РЕТМТХ, l 
Jc cech e he he hh he e khe ke hk hee hehe | 


i 
PRINT*| ' 
PRINT*, I * жі 
PRINT*, ! ж 1. MULTIPLE COLUMN FORMATION *! 
PRINT*, 1 * *! 
PRINT*, j * 2. SINGLE COLUMN FORMATION ж! 
PRINTS, 1 * ж! 
РВТМ 1 kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkki 
SEIT ж) UNITE 
PRINT*, κα NOTE : FOR YOUR REFERENCE, CURRENT ANSWER IS', 
T 

а. ( S TIR -EO. 1 ΟΡ. UNITF ο ГЕЦ 

ο 

ELSE IF ( UNITF .NE. 1 .OR. UNITF .NE. 2 ) THEN 
PRINT*, ! 
PRINT*, '*** ERROR : ENTER THE NUMBER 1 OR 2 *Х*х*! 
PRINT*, ' 

GO TO 17 

ο 55 
WRITE(6,603) 
FO BMAD el ω. 
PRINT*, ! 
PRINTT, | DO του WANT TO USE SPECIAL RULE FOR BOUNDARY?! 
PRINT*, | YES ( 1 ) OR NO ( 2 )! 
READ(5, ж) NSPCL 
PRINT*, !** RM : FOR YOUR REFERENCE, CURRENT ANSWER 15', 

L 

WRITE(6,604) 
FORMAT, ШРІ) 
PRINT*, | WHICH OUTPUT DO YOU WISH TO PRINT? ( 1 -5 )'! 
PRINTS. i Z e de de e de dee v ee ! 
PRINT*, | ж 1. GENERAL INFORMATION FOR EACH ARC ж! 
PRINT*; | ж 2. FLOW RATE & TRAVERSAL TIME ты 
PRINT*, ' * 3. INPUT DATA FOR CHECK κ! 
PRINT*, ж 4. COORDINATES OF BOUNDARY LINE ж! 
РЕТМТХ, | ж 5. NETWORK DATA FOR NETWORK PROGRAM ж! 
PRINT*, 1 kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk i 
READ(5, ж) РКІМТ 
PRINT*, '** NOTE : FOR YOUR REFERENCE, CURRENT ANSWER IS', 


PRINT 
IF ( PRINT .EQ. 1 .OR. PRINT .EQ. 2 .OR. PRINT .EQ. 3 
OR. PRINT .EQ. 4 .OR. PRINT .EQ. 5 ) THEN 
GO TO 19 
ELSE IF ( PRINT .NE. 1 .OR. PRINT .NE. 2 .OR. PRINT .NE. 3 
те .OR. PRINT .NE. 4 .OR. PRINT .NE. 5 ) THEN 


PRINT*, '*** ERROR : ENTER THE NUMBER 1 THRU 5 ***! 
PRINT*, | | 

GO TO t8 

END IF 
THRESHOLD FOR BARRIER DUE TO SLOPE OFF ROAD 

THRES1 0.3 


ТНБКЕ52м-енбке 


x 2. SET THE INITIAL CONDITIONS. 


А TOTAL NUMBER OF ARCS IN NETWORK. 
А TOTAL NUMBER OF NODES IN NETWORK. 


NODE 
CALL INIT(NX, NCODEX,XL,YL,XR,YR,XA,YA,SLPCH,HI,NO,DHTI,CASE, 
DISTA, SLOPEL, ЛЕР ХБ” ‚НТ2, DHT2, рІЗТБ. SLOPE2, CODE) 


NARC 


* 3. READ THE DATA FROM THE DATA FILE. (DATA 1, DATA 2, DATA 3, DATA 4) 
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κ DATA 1 ; ALTITUDE AND CODE 
κ DATA 2 ; NODE CHARACTERISTICS 
κ DATA 3 ; ARC CHARACTERISTICS 
κ DATA ; SPEED FOR MOVEMENT 
CALL INPUT1(NARC,NODE,ALTD,NCODE,XNODE,YNODE, 
ж UNITE, CASÉ, ТНВЕ51, THRES2, SP ,PRINT) 
NCOUNT = 1 
ІҮ = 100 
DO 10 ТВ = 1, МАКС 
| CALL INPUT2 (XNODE, ο ος ο ШІ ІХ2, ГҮ2 МОРЕА МОРЕВЕ 
~ БАСІР ШІП ÜTMY) 


* 





* 4. DETERMINE THE ANGLE BETWEEN ARC AND X-AXIS. 





БЕНТАЯ = TX2 Tr] 
DELTAY = IY2 - ТҮ1 
IF ( DELTAX .NE. 0.) THEN 
XM =M DELITAY / DELTAX 
ELSE ы ( DELTAX Ж О.) ТНЕМ 


= 2999998 
END IF 
RADIAN = Soe ХМ ) 
DISE ο οτι DELTAX 5“; 2 + DELTAY ** 2) 
ANGLE] = ( RADIAN ^ 180. Me e τ 141592 


F ANGLE! „СОЕ. О T 
ANGLE = ANGLE1 
ЕНСЕ τα TANGLET LT. 0.) THEN 
ANGLE = 180. + АМСІЕ1 


END IF 
ж 
* 5. COMPUTE THE DISTANCE BETWEEN EACH PAIR OF POINTS ALONG ARC. 
ж 
ΠΤΙ τ; ο ο Ισ , DIST, <M,XA,YA,DISTA,K) 
τ 6. COMPUTE AVERAGE SLOPE ALONG THE ARC. 
LX1 - ІХІ / 100 
үй = νι ος 
LX2 = IX2 / 100 
Буг = ТҮЗЕ /Й 100 
ALT1 = ΠΠ РАЈЕ Dd 
ALT2 = Αυτο κο. LY 
AVGSLP = ( ALT2 - ALTI ) / DIST 


А 7. COMPUTE SLOPES BETWEEN EACH PAIR OF POINTS ALONG THE ARC. 


CALL SLOPI(ANGLE ,ALT1L,ALT2,IX1,IY1,1IY2,1X2,DISTA,XM,XA,YA,ALTD, 
Ж aT DBT SLOPE Nay ACTD AD) 


* 8. AS A MEASURE OF TERRAIN " STEEPNESS " ALONG THE ARC, COMPUTE 
А TOTAL AND AVERAGE SLOPE CHANGE. 
ASLPCH = 0. 
TSLPCH = 0. 
NB = y: τ 
το 11 = ΝΒ 
SLPCH(ÍD) = ABS( E - SLOPE1(ID)) 
TSLPCH = TSLPCH + SLPCH(ID) 
ү CONTINUE 
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ASLP 


9. COMPUTE SLOPE 


CALL SLOP2 


CH = TSLECH / NE 


S FROM ARC TO POINT OFF ARC, AND MIN DISTANCE. 


* 10. DETERMINE COORDINATES FOR THE MIN DISTANCE.( LEFT AND RIGHT ) 


CALL MIND(K,XM,YM,XA,YA,LMIN,R8MIN,XL,ZL,XR,TR) 


* 11. COMPUTE THE RATE OF FLOW ALONG THE ARC. 


CALL FLOW(S 


* 12. ος ασ THE 


P,WIDTH,ARCTP, CASE ,FR,MAXNO) 


MIN TRAVERSAL TIME ALONG THE ARC. 


CALL TIME(ARCTP,SP,UNITE ,CASE, DIST ,WIDTH,ATIME) 


Ж 13. ӘРКІМ ІНЕРСЕ 


CALL РКЕ 
NCOUNT = NC 

10 CONTINUE 
STOR 


END 
ΧΧΛΧΧΧΛΑΧΧΧΧΑΧΑΧΧΧΧΧ 


* A. SUBROUTINE FO 


Хи 


„ SUBROUTINE 


INTEGER 

DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 


SULT OF CALCULATION. 


OUNT * 1 


wk xx 


R SETING THE INITIAL CONDITIONS. 


1 < Z Z Z Z Z Zç Zç Z Z Zç Zç Zç Zç Zç Zç c Zç Zç Xç Zç Zç Zç Zç Z< Z< Zç Zç Zç Zç Zç Zç Zc Zç Zç Zc Zç Zç Z Zç Zç ç k k k t k kk k k k 


INIT(NX,NCODEX, XL, YL,XR, YR, XA) Ya, SLPCH HT NO, Dir Case 
DISTA, SLOPE1, Авив ,HT2, DHT2, DISIE? SLOPE2 , CODE) 


CASE,CODE 
XL(10) ТОЕ M 78(10), NX(10),NCODEX(10) 
DISTA(90),XA(90),YA(9 δ) 

НТ(90), DHT (90) , SL о 

DISTB(100,100 xb 190), ,B(100,100) 
HT2(100,100),DHT2(100,100) SLOPE2(100,100) 
NO(90), SLPCH(90 ) 

CODE(0:7,0:7) 


= INITIALIZATION S Mm d 


DO ІС E 
XL(L 
YL(L 
AR(L 
TR(L 
NX(L 
NCOD 

10 CONTINUE 
DO. TIMES 


SLPCH z 
A (L 


DISI 


, 


- 


) 


) 


ООООО 
ο 2 ο H+ 


| 
о 


EX(L) 


r^ rptu ii 


Pri nud 


{| 
ООООО «0 


| | • 
© 


| 
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SLOPE1(LA) = 
II CONTINUE 
ΡΟ 12 LB = 1l, 100 
BO -I2 LE = 1 


ы 
Hj 
N 
үтчет- 
r< 
UJ - 
4-2-7 
(а 
Νο 
H H H HI H Hl 
ΟΟΟΟΟΟ 


ОБОРЕ ТВ не 
12 CONTINUE 


x SET CODE OF BOUNDARY FOR VEHICLE UNIT. 
τ CASE ps 1) THEN 


00 13 Ј = 
CODE(I,J 
4 


DISTB T ds 


~ 
| 


(а 

O 

о 

[т] 

fs 

-3 

О О 

uw i ` ~ 
PB o 


13 CONTINUE 


С Dm 
х SET CODE OF BOUNDARY FOR DISHOUNTED TROOPS. 
5р Я E μμ 2) ΤΗΕΝ 


DO 14 J 
CODE CE, 3) z 
14 CONTINUE 

DO l- sk 


0, 
CODE( IA,L 

CODE(IA,3 

GODECIA,5 

15 СОМТІМУЕ 


, 


ES 


uut 
ооо 


END : 
kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


* B. SUBROUTINE FOR READING THE DATA FILES( DATA 1, DATA 2, DATA 4 ) 
kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


„SUBROUTINE INPUT1(NARC,NODE,ALTD,NCODE, XNODE , YNODE. 
UNITF ,CASE, THRES1, THRES2, SP , PRINT) 


INTEGER ALTD,CASE,PRINT,UNITF 

CHARACTER  *1 BLANK 

DIMENSION ALTD(100,100) ,NCODE(100,100 

DIMENSION B ,NCODEX(10), XNODE(90), )wopE (90) 
DIMENSION 5Р(7,3,2 


* READ THE ALTITUDE AND CODE BY TWO DIMENSIONAL ARRAY. 
DO 20 IN - 1,99 | 


Т =1 
28-81-11 
pO 21 ІТ = 1,10 
READ(1,100) ( El) NCODEX(I1), I121,10 ) 
100 FORMAT(I4,12,9(15,12)] 
τυ ΤΟΞΟ 
ALTD(I,J) = NX(IS) 
NCODE(I,J) = NCODEX(IS) 
Т = + 1 
22 CONTINUE 
од CONTINUE 
READ(1, 101) BLANK 
101 FORMAT (A1) 
20 CONTINUE 
* READ THE NODE CHARACTERISTICS FOR SINGLE ARC ATTRIBUTES. 
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DO 23 IA = 1,NODE 
READ(2,11 10) ,ENODE (1A) ,YNODE (1A) 
110 FORMAT( 4X,2F4.0) 
23 CONTINUE 
DO 26 IM s 1,7 
DO 26 JM s 1,3 
DO 26 КМ = 1,2 
* READ THE SPEED FOR MOVEMENT BY THREE DIMENSIONAL ARRAY. 
READ(4,130) SP(IM,JM,KM) 
130 FORMAT(F5.0) 
26 CONTINUE 
IF ( PRINT 1 ТНЕМ 
НЕІТЕ (6, 7167656, UNITF,THRES1, , THRES 2 
716 FORMAT ( ' 4X, UNIT TYPE 4,/,5X,' UNIT FORMATION', 
* ig! m ,/,5X,'POS. THRESHOLD a F7.3, 
κ / 5X, NEG. ‘THRESHOLD ΠΠ) 
ELSE ΙΕ ( РВІМТ EQ: 2 THEN 
WRITE (20,7 3) CASE, UNITF,THRES1, , THRES2 
769 FORMAT ( ! 4X, UNIT ТҮРЕ ,14,/,5X,'UNIT ЕОВМАТТОМ', 
ж | ,14 ,/,5Х%,'Р05. THRESHOLD :!, Е7.З3, 
ж /,5Х'. 'МЕб. ‘THRESHOLD ,F7.3 ) 
ELSE IF ( PRINT 1199 3 ) THEN 
WRITE(6 
710 FORMAT(! “л 71('-'),//,4X,' INPUT DATA !, 
κ M ALTÌTUDÉ AND CHARACTERISTICS )',/) 
DO 24 NA = 1,2 
DO 24 NB = 1,100 
WRITE(6,711) NK,NB,NA, ALTD(NB,NA) ,NCODE(NB,NA) 
a FORMAT(! 15,4 (1, 1321570: ШЕГЕН 
7 NK = NK + 1 | 
24 CONTINUE 
WRITE(6,712) 
п TES ,3X,//,4X,' NODES ALONG ARC ', 
WRITE(6,713) (XNODE(IA), YNODE (1A), A=1,,NODE) 
713 ҒОРМАТ(! !,6Х, 14,3Х,2Е7.0) 
ELSE IF ( РВІМТ .Е0. 5 ) ТНЕМ 
МЕТТЕ (6,71 
714 FORMAT (' ' 6X,'HEAD TAIL',3X,'TIME FLOW RATE DIST WIDTH SPEED!) 
I 
RETURN 


END 
ACA hh CH C ACA 


k C. SUBROUTINE FOR READING THE DATA FILE ( DATA 3 ) 
kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


SUBROUTINE INPUT2(XNODE , YNODE ТҮ,ТХ1,ТҮ1, IX2,IY2,NODEA,NODEB, 
ARCTP, UTMX , UTMY ) 


INTEGER ARCTP , UTMX , UTMY 
DIMENSION  XNODE(90),YNODE(90) 


< READ THE ARC CHARACTERISTICS FOR SINGLE ARC ATTRIBUTES. 
READ(3, 1 a): NODEB , ARCTP 
120 РОВМАТ (13,1 
МХ] E RODE NODEA | 
NY1 = YNODE( МОРЕА 
UTMX= NX1 / 100 
UTMY- NY1 / 100 
UTMX= UTMX * 10 
UTMY= UTMY * 10 
МХ1 X МА IY ) * 100 
MZ 1 MOD( NY1,IY х 100 
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NX2 ee NODEB 

NY2 - YNODE( NODEB 

ΜΧ2 Wes. NX2,1Y * 100 
MOD(. NY2,IY ) * 100 


MY2 
IF ( MY2 .GT. MY1) THEN 
MX1 


TAL 
Izi МҮ1 
152 ΜΧΖ 
IY ΠΕ | 
% CHANGE THE ORDER OF NODE FOR USING ANGLES OF INCLINATION. 
peor ( MY2 217. MYl)XTHBN 
ІХ] МХ2 
T МҮ2 
T М 
MYI 
NODEA 
NODEB 
NODE2 
το, 
ECSETIE } Mize 5 ean THEN 
TE MX2 .GE. MX1) THEN 
- МХ1 
= МҮ] 
1Х2 - МХ2 
а = Мег 
БЕЗЕР M2 «Lim MXP THEN 


224 
О 
О 
mi 
"n 
H H H H 


DTI 


14] Ic 
ТҮЙ МҮ2 
T МХІ 
τα MYi 
NODEA 
NODEB 
ΝΟΡΕ 2 
МОРЕ] 


= 
О 
Oo 
4 pd 
D? 
H H B Hl 


END IF 
RETURN 


END 
T Z Z Z Zc Z Zc 27 Z Zc Zc 7 7k 7 7 Zc 7k 7k Z Zc Ze Ze Ya 7k 7 Ze yk 7k Zc Z Zc Zc ye 7 7 Zc Z Z Zc Zc Zc Zc Z e 7 Z Zc Zc xk dx xk Ye 7k $k Zc 7 Zk Z Ze Z 7k e Zc Ze Zç Ze 


* D. FUNCTION FOR CALCULATING GRID COORDINATES( X ). 
TX X Yo Zç 7 yc 7<ç Zc 7 Zç Zc Tç Zc Zc Zç Zc Zç Tc Zc Zc Tc Zc Zc Tc Ye Y Zc Z< Zc Zc Zc Zc Zc Z d$k dede Zc Zc 7 Zc Zc Zç 7ç 7 Zç ç Zç Yc ç 7ç Zç 7 Zc Zc Z Z Zc Z< Tç Zç Z Z =< 
FUNCTION XCORD(XM,IA,IY,IX) 
IF ( XM .МЕ. 0. ^ THEN 
XCORD = (( 1. ХИ ) ж к TA - IY ))+ IX 
ELSE IE IE. 4 ΧΜ E: ο TH 
CO 


ІХ 
ЕМр IF. 
RETURN 


END 
xxi xk ЈОЈ eek 


x &. FUNCTION FOR CALCULATING GRID COORDINATES ( m. 
CCCo eee ee e ЛО c k k k k k kk ok kk k k ЖЖ 
FUNCTION YCORD(YM,XXB,YC,XC) 


οσο = ЧУМ ΤΕΒΕ = ACs) \ ҮС 
RETURN 


END 
kkkkkkkřžkřkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


* F, SUBROUTINE FOR COMPUTING THE DISTANCE BETWEEN EACH PAIR 
ж OF POINTS ALONG ARC. 


ΧχΧΧΧΧΧΧΧΧΧΧΑΧΧΑΧΧΧΧΧΧΧΧΧΧΑΧΧΧΧΧΑΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΑΧΧΑΧΧΧΧΧΧΧΧΧΧ κ 
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SUBROUTINE DST(IY1,IX1,IY2,IX2,ANGLE,DIST,XM,XA,YA, 
* DISTA 


‚К) 
DIMENSION DISTA(90),XA(90),YA(90) 


IF ( ANGLE .GE. 45. .AND. ANGLE .LE. 15, ) ТНЕМ 
МА - ((ІҰ2 - 191 + 49 ) / 100) + 1 
ELSE IF ANGLE .GT. 135. .OR- ANGLE .LT. «5; ΤΠ 
NA4 = I = АЊА 49 ) / 100) + 1 
NA135 = LAL τι τ» 
a ANGLE .LT. 45.) NA = N 
IF( ANGLE .GT. 135.) NA = ΝΑΙ55 
END IF 


K = 
ЧА(К) = 191 
ХА(К) т IXI 
~ ANGLE OF TNCLINATION IS BETWEEN 45 AND 135 DEGREES. 
IF ( ANGLE «68. 45. .AND. ANGLE .КЕМІЗ5. ) ІНЕП 


DO 30 IA = 0, IY2-100, 100 
DENTS SERV 
IF (A .GT. 0.) THEN 
К = К+1 
Tye =e 
фер 550 
ШЫ = REAL( ІА ) 
ХА(К) = XCORD( XM,IA,IY,IX ) 
END IF 
30 CONTINUE 
x ANGLE OF INCLINATION IS GREATER THAN 135 DEGREES. 
ELSE IF ( ANGLE .GT. 135. ) THEN 
E | IXl + 100 ) / 100 
12 = ( 12 / 100) + E 
Ρο 31 IA = 11, nos 
A(R) = IA * 100. 
MONS TUM 
КС = ТА 
ҮМ = ХМ 
ХХВ - ХА(К) 
YA(K) = YCORD( ҮМ, ХХВ, YC, XC ) 
31 CONTINUE 
* ANGLE OF INCLINATION IS LESS THAN 45 DEGREES. 
ELSE IF ( ANGLE .LT. 45. ) THEN 
iL e IXl + 100 ) / 100 
ӨЛЕ бо ако) ~ 1 
DO αρ PA τα O 
= K + T 
ХАК) = ТОТИ 
σοι τν 
πο ia 
УМ = ХМ 
ХХВ = ХА(К) 
YA(K) = YCORD( YM, XXB, YC, XC ) 
32 CONTINUE | | 
END IF 
ж COMPUTE THE DISTANCE BETWEEN БАСН PAIR OF POINTS. 
aw CSS] 
LEON 
R τη 
D cn 
IF( XM .NE. 0. ) THEN 
YA S = 1/2 
ХА(К) - XCORD( ХМ, IA, IY, IX ) 
ELSE IF( XM .ЕО. 0. ) THEN 
YA Е zs 
ХА(К) = IX2 
END IF 
ртт б 
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DO 33 ІВ = 1, МА 
IF( IB .LT 


NA _) ТНЕМ 
DISTA(IB) = SQRT( D nep - D хх2 
ж ш YA(IB*1) - YA(IB)) **2 ) 
ELSE τος ІВ -EQ. A ) THEN 
ТА(ТВ) = DIST - DISTI 
ND τι 
DIST1 = DIST1 + DISTA(IB) 
33 CONTINUE 


RETURN 
END 
ek sk heh khe he e e e i e e e he he e eh he e he he he e he heh e e hehe he he e he e he he hh e he e he he heh he he hh hehehe hehe 


ж G. SUBROUTINE FOR COMPUTING SLOPES BETWEEN EACH PAIR OF POINTS 
* ALONG THE ARC. 


3 kk e ve e ve v e e ve e e he e e ehe v e e he hee e the he e zh he e ehe e e hee e he e hehe e e cx 


„SUBROUTINE SKLOPI (ANGLE ΕΤΗ ΑΙ το ΙΕ IY1 IY2 IX2 DISTA, RM KA,YA, 
ALTD, HT, DAD, SLOPE1, NA, ACTD ,AD) 


INTEGER ALTD 


DIMENSION ALTD(100,100) DISTA(90) XA(90), YA(90) 
DIMENSION HT(90) ,DHT(90) ,SLOPE1(90) 
DIMENSION ACTD(90) 
AD = 0. 
RY = 100. 
HT(1) = ALT1 
IF (ANGLE .GE. 45. .AND. ANGLE „ТЕ. ν THEN 
МА ((I¥2 = Tvl + 49) / 1 coe 
DO 40 Ic = 2, МА 
ЈХ1 = INT ( Ka (TC) / 100 
JY1 = INT( YA(IC) / 100 
R = MOD (HA (TC) . RY 
HIA - (В/100.) * ( ALTD(JX1,J3¥1) - ALTD (JX1+1, ЈУ1)) 
CIC LT. NA ) THEN 


END 
DHT. ICSI = HT(IC )- HT(IC-1) 
SLOPE1(IC-1) = DR fc^ 1) 4 E -1) 
ACTD(IC- ) = DISTA(IC-1) / COS(SLOPE1(IC-1)) 
AD = AD * ACTD(IC-1) 
40 CONTINUE 


ELSE IEI ANGLE TGTS AROR: ANGLE .LT. a ) THEN 
= 





МА45 - 114 - IX1+ 49 | / 100) 
NA135 = ((IX1 - IX2+ 49 / 100) + 5 
ANGLE .LT. а МА = МА45 
IF( ANGLE .GT. 135.) NA = NA135 
DO 41 IC = 2, NA 
J2 = x e / 100 
JY2 - τς % 100 
R T GE Έτος ҢҮ) 
ΗΤΑ = (R/100.) ( ALTD(JX2, JY2) - ALTD(JX2,JY2+1) ) 


1 ICN ET. ΠΕΙ 

HT(IC) » ALTD(JX2,JY2) - НТА 
ELSE IF( ТОЙ БО. NÀ ) THEN 

НТ(ІС) = ALT2 
END 


DHT(IC-1) HT(IC)- HT(IC-1) 
SLOPE1(IC-1) DHT(IC-1) ISTA (16-1) 
ACTD(IC-1) DISTA(IC-1) / COS(SLOPE1(IC-1)) 
AD = AD + ACTD(IC-1) 
41 CONTINUE 
END IF 
RETURN 


END 5 
kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 
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* H. SUBROUTINE FOR COMPUTING SLOPES FROM ARC TO POINTS OFF ARC. 
kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


72 


SUBROUTINE 
INTEGER 
REAL 


DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 


SLOP2 
ELD ,ARCIP,CASE,CODE,PRINT,UIMZ,UTMY 


ALTD(100,100),NCODE(100,100 NX (10) ,NCODEX(10) 
BUT. „“МОРЕ (90) ,СОРЕ (0:7,0:7) 

DISTA(90) ,XA(90) , YA(90 

HT(90), DHT (30) , SLOPE (90 

DISTB( 100 100 m oe , YB( 100,100) 
HT2(100, 100) , DET2 100,100) ,SLOPE2(100,100) 


NO(90) ,SLPCH(90) 
XL(10),YL(10),XR(10), I EM, 
SP ο... ,MAXNO(3,3) 
C IME(3 

ROADW( 


COMMON / INTS / NOTER IRI ISl Ime le ES HAS UE 


DEA,NODEB,J1,J2,J3,J4, PRINT;UTIMX,UTIITOS 
CASE,ARCTP,FR,AD,UNITF,NSPCL 


COMMON / INTA / NCODÉ , NO, ALTD, CODE , MAXNO 


COMMON / SR 
x 


EAL / XM, ANGLE, DIST, TSLPCH, ASLPCH,ALTl,ALT2,AVGSLP, 
LHIN, RMIN, AL, Y ZR, УК, SII WIDTH, THRES1, THRES2 


д COMMON / AREAL / A YA, Н; DHT , DISTA ,DISTB, SLPCH, SLOPE1, SLOPE2, 


ВОВ; НТ2, DHT2, sP. АСТР, АТТМЕ, ROADW 


LMIN = 10000. 
RMIN = 10000. 
= 100 


RY 


ANGLE OF INCL 
IF ( ANGLE 


KD 
IH 
Τη 
ХВ(І 
УС 
ХС 
ХХВ 


ТК XB(T 


END 


YM 1. / ХМ 
ELSE IF ( XM EQ. 0.) THEN 
YM - 9999999. 
END IF 
INATION IS BETWEEN 45 AND 135 DEGREES. 


„СЕ. 45. .AND. ANGLE .LE. 135. ) THEN 
COMPUTE THE S 
. DO 70 IE 


B EQUI THE LEFT SIDE. 


, 


= 1 

= (INT(XA(IE) / 100) * 100)- 100 

E,KD) = ІН 

= ЧА(ТЕ 

= ХА(ТЕ 

= XB(IE,KD) 

E,KD) = ZCORD( nt, XXB, YC, XC) 

= ЕЗІ Жр, у 

= YB(IE,KD) / 100 

= ТВ ТЕ КО) Ба RY)) / 100. 
ze RA СЕ 06 A ) J4 = J2 + 1 
ΤΕ. RA LIT EY: 5 24 - 22 

ELSE IF ( XM .LT. ) THEN 


ТЕ( ВА СЕ. О : b IHEN 
J4 = 22 
ELSE E RB. .LT.. 05 ) THEN 
24 - 22-1 
END IF 
END IF 
E; КИ .GT. 10000. .OR. YB(IE,KD) .GE. 10000. ) THEN 
D 


ТЕ 
DISTB(IE,KD) = SQRT(( XB(IE,KD) - XA(IE)) **2 


70 


78 


ж 


ж 


74 


ТЕ ( 


ELSE 


+ ( YB(IE,KD) - YA Bn х2) 
ВА 4 (аіто(21,74) - - ÁLTD(J1,J2*1 
D(Jl B 


NCA - 7 
PESE TF. KD POr le) THEN 
2А т XB (TE, E / 10 
= YB(IE,KD-1) / 10 
МСА = МСОРЕ( ЈА,ЈВ) 


END IF 
NCB = NCODE(J1,J4) 
CODE(NCA,NCB) .50. 1 .AND. SLOPE2(IE,KD) .LE. ТНВЕ51 
ТШ „SLOPE2 (IE, KD) .GE. THRES2 .AND. KD .LE. 10) THEN 
= I 


ТЕ (GORE NGA НСВ | ПЕ. 1 .OR. SLOPE2(IE,KD) .GT. THRES1 
„ОВ. SLOPE2(IE, KD) јр THRES2 .OR. KD .GT. 10 ) THEN 
ІН - ІНІ % 200 

IF(DISTB(IE,KD) .LT. LMIN) LMIN = DISTB(IE,KD) 


APPLY SPECIAL RULE (WIDTH DETERMINATION RULE 2) FOR BOUNDARY. 


IF( NSPCL .EQ. 1 ) THEN 
IF( KD .EQ. 1 ) THEN 
ІЛІМ = 0. 


ENDE IF 
END IF 
GO TO 71 
END IF 
COMPUTE THE SLOPE ON THE RIGHT SIDE. 
KD = KD + 1 
ХЕ) ТЕ, AD) = ІН 
= BUE 
ae = XA(IE 
ХХВ = XB(IE,K 
YB(IE,KD) = Z ORDE. a ΠΕ, τς, ΠΕΠ) 
ЈАК S SIP TE RD) 7 1 
J2 = YB(IE,KD) / on 
ЕВ = (MOD (YB(IE,KD), RY)) i 100 
D2 = ABS( YA(IE) - YB(IE,KD) ) 
IP ΧΜ Е. 0. ) ТНЕМ 
a RB .GE. 0.5 ) J4 = J2 + 1 
TP το Ὁ Ὁὓ,ς „AND. D2 .СІ. 50.) 74 - 72 
ELSE ТЕ πω. ο. ) THEN 
= RB .GE. 0.5 ) ЈА = Ј2 + 1 
τ МЕН 115 0.5 ЈА = J2 


END 
Ἢ .GT. 10000. .OR. YB(IE,KD) .GE. 10000. ) THEN 


END 


IF(DISTB(iE, KD-1) .LT. RMIN) THEN 
RMIN DISTB(IE, KD-1) 





END IF 
GO T0877 
DISTB(IE,KD) = SORT ( ( EB (IE. RD] - А (ТЕ | хж2 
YB(IE,KD) - ҮА(ТЕ)) х%2) 
HTB = RB * (ALTD(J1,J2) - ALTD(J1,J2+1) 
HT2(IE,KD) - ALTD(J1,J2) - HTB 
DHT2(IE,KD) = HT2(IE,KD) - HT(IE) 
IF( DISTB(IE,KD) „EQ. 0. GOR TO 
SLOPE2(IE, хр) = DHTZ(IE, Ὃν / DISTB(IE,KD) 
ТЕЙ Т EQ, IH1+200 ) THE 
ELSE TF ( KD. ӨШ І ) ТНЕМ 
= . KD- 1 / 100 
JB = YB(IE,KD-1) / 100 
NCA = NCODE( JA,JB 


END 
NCB = NCODE(J1,J4) 
71 


76 ТЕ ( СОРЕ(МСА,МСВ) .Е0. 1 .AND. SLOPE2(IE,KD) .LE. THRES1 


* .AND. SLOPEZ(IE, KD) „СЕ. THRES2 .AND. KD .LE. 20) THEN 
ІН = IH + 10 
KD = KD + 1 
GO TO 74 


ELSE IF (CODE(NCA,NCB) .NE. 1 .OR. SLOPE2(IE,KD) .GT. THRES1 
‚ОК. SLOPE2(IE RD) БИ THRES2 .OR. KD .GT. 20) THEN 
IF (DISTB(IE, KD) LI. AMIN) THEN 
RMIN DISTB(IE KD 


END IF 
Е APPLY SPECIAL RULE(WIDTH DETERMINATION RULE 2) FOR BOUNDARY. 
IF( NSPCL .50. M THEN 


IF( KD .EQ. 2 ) THEN 
RMIN = 0. 


END IF 
END IF 
KD = KD + 1 
END IF 
77 NO(IE) = KD - 1 
70 CONTINUE 
* ANGLE OF INCLINATION IS GREATER THAN 135 OR LESS THAN 45 DEGREES. 


ELSE IF ( ANGLE .GT. 135. .OR. ANGLE .LT. 45.) THEN 


^ COMPUTESTHE SLOPE ΟΡ ΕΕ ШЕК Ж ер» 
DO S90 ЕЕ IK 
KD = I 
ШЫ ANGLE .GT. ο ΤΗ 
LE In 


62 YB(IE,KD) = IH 


uc do А 100) - 100 
ІМТ(УА(ТЕ)/100) * 100) + 100 


E, KD) = XCORD( YM, ІА, ІҮ, ІХ) 
( XB =) THEN 


ТЕ 
ae 2) 1 ΤΠ 
YB(IE,KD) / 1 
(προ. ΕΤΕ RD), ΠΟ ТОСЕ 
IF( RA EL 
J3 = E τὶ 
ELSE IF( RA T. 0.5 ) THEN 
J3g- JT 


END IF 
IF ( RB QE Tc 10000. .OR. YB(IE,KD) σε. 10000. ) THEN 


END IF Ë 
DISTB(IE,KD) = SQ xB E Ts: 20) т СЕ 
+ ҮВ КО) - ҮА(ТЕ 
НТВ = БА Х а> - ALTD(J1+1, J2 
HTZ(IE,KD) = ALTD( J1 J20- IIB 
DHT2(IE, KD) = HTZ(IE, FOU, - 1) 
ТЕП DISTB(IE, KD) τη P) GO TO 65 
SLOPE2 (TEF Кр) = DHT ГТЕ? RD) / DISTB(IE,KD) 
ТЕ ( KD T Е ) THEN 
C = 


ELSE τες KDR CTED THEN 

ЈА = ЕЕ: „Дре i / 100 
= YB(IE,KD-1) / 100 

NCA = NCODE ( JA,JB) 


_END ТЕ 
NCB = NCODE(J1,J4) 
65 IF( ANGLE F OTR) THEN 
ЕНИ CODE(NCA, NCB) «ΕΟ. 1 .AND. SLOPE2(IE,KD) .LE. THRES1 
"D SL0PE2 (E Κδ) „СЕ. THRES2 .AND. KD .LE. 10) THEN 


κκ 2) 





| х*2 
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AB EE no s 1 
ELSE τ “(CODE (Nea, ПРЕ. ‚МЕ. 1 .OR. SLOPE2(IE,KD) .GT. THRES1 
e εις. Б. КЕЙ .LI. THRES2 .OR. KD .GT. 10) THEN 
= + 


T 20 
IF(DISTB(IE,KD) .LT. LMIN) THEN 
ue = DISTB(IE, KD ) 


END 
х APPLY SPECIAL RULE (WIDTH DETERMINATION RULE 2) FOR BOUNDARY. 
IF( NSPCL EQ. 1 ) THEN 
IF( K D +E. 1 ως THEN | 


END IF 
END ΤΕ 
GO TO 61 


END IF 
ELSE IF( ANGLE .LT. 45.) THEN 
IF ( CODE(NCA, NCB) ON АМЫ SEDOPEZ(IE,KD) .LE. THRESI 
.AND. SLOPEZ(IE KD) .GE. THRES2 .AND. KD .LE. 10) THEN 


pee Le + 
KD = KD + 1 
GO TO 62 


ELSE SECO ENED Мав) Γρ OR. SLOPE2(IE KD) СТ. THRES1 
. SLOPE2(IE , KD) .LT. THRES2 .OR. KD .GT. 10) THEN 
E ІНІ - 200 
ΠΕ KD) ELIT.: LMIN) THEN 
LM = DISTB(IE, KD) 


END ТЕ 
* APPLY SPECIAL БЕРЕ, НІНЕН Орал ыы RULE 2) FOR BOUNDARY. 
IF( NSP dub. THEN 
TF( RD Eo. 1 Ж ШОҒЫ 


END IF 
END IF 
GO. TOwGI 


* COMPUTE THE SLOPE ON THE RIGHT SIDE. 
61 KD = KD + 1 


À m 
ХА(ТЕ 
ШЕКТЕ, КІ) 
KD) = ODE ο. IA, Pe 1X) 
КЕЕ, ,KD) / 1 
YB IE КОИ 100 
ME (4B (TE, 0:57). ВЕ / 100. 
LEO RB .GE 5 THE 

J3 E 1 
а ББ 114. 0.5 ) ТНЕМ 


> 
UJ 
~ 
H 
II l lc H H I 


END 
IF(XB(IE, KD) ТОТ. ТОООО ОВ. YB(IE, E .GE. 10000. ) THEN 
IF(DISTB(IE, KD-1) es RMIN) THEN 
RMIN = DISTB(IE,KD -1) 


END IF 
59 Т0 67 
ЕМІ) 
DISTB(IE, KD) = ФЕ м 2 - E ) хх2 
ж YB(IE,KD) - YA(IE)) **2) 
НТВ = RB * (ALTD(J1,J2) - ALTD(J1*1,J2) 


HT2(IE,KD) = ALTD(Ji,J2) - НТВ 
DHTZ2(IE,KD) = HT2(1E, KD) - HT (TE) 
IF( DISTB(IE,KD) EQ, 0 GO TO 66 
510РЕ2 (ТЕ, KD) = gne (IE,KD) / IDISTB(IE,K D) 
IF( ANGLE .GT. 135 .AND. IH .EQ. IH1+200) THEN 


NCA = 7 
ВИЗЕ art ЕТЕ LT. 45 .AND. IH .EQ. IH1-200) THEN 


73 


66 


67 
60 


ELSE IF( KD .GT. 1 ) THEN 
JA = ΕΕ} ,KD- E / 100 
JB = IE,KD-1) / 100 

NCA = NCODE( JA,JB) 


END 
NCB = NCODE(J1, J4) 
IF ( CODE(NCA,NCB) 5), .AND. SLOPE2(IE,KD) .LE. THRES1 


.AND. SLOPE2(IE .GE. THRES2 .AND. KD .LE. 20) THEN 
IF( ANGLE .GT. 135.) THEN 
IH = IH : 100 
END IF 
IF( ANGLE .LT. 45. ) THEN 
IH = IH - 10 
END IF 
IF( IH ΤΕ. 0 ) GO TO 60 
KD D. KD, 
ELSE IF) CODE (NCÀ, NCB) 1 .OR. SLOPE2(IE,KD) .GT. THRES1 


ОК. SLOPE2(IE, KD) NI. THRES2 aves KD .GT. 20) THEN 
TF (DISTB(IE, KD) υπ. ЕМІМ) Т 
КМІМ - DISTB(IE, KD) 
END IF 


APPLY SPECIAL RULE (WIDTH DETERMINATION RULE 2) FOR BOUNDARY. 
IF( NSPCL EQ. THEN 
F( K B o P 2 Ji THEN 
IN 


END IF 
END IF 
KD = KD + 1 
END IF 
NO(IE) = KD - 1 
CONTINUE 


END IF 
DETERMINE THE TOTAL ARC WIDTH. 
ASSIGN THE WIDTH OF ROAD ITSELF. 
ROADW 18. 
ROADW 
ROADW 
ROADW 
ROADW 
ROADW 
ROADW 
IF (ARCTP .NE. 7) THEN 
WIDTH LMIN, + RMIN + ROADW(ARCTP) 
ELSE IF (АВСТР. .Е0. 7) ТНЕМ 
WIDTH = ROADW(ARCTP) 
END IF 


RETURN 


SOUP WN EH 


6 
4. 
4. 
2. 
5 
4. 


END 
kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


* I. SUBROUTINE FOR COMPUTING COORDINATES FOR THE MINIMUM DISTANCE. 
kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


SUBROUTINE MIND(K,XM,YM,XA,YA,LMIN,RMIN,XL,YL,XR, YR) 


REAL LMIN 
DIMENSION SMS и 
Ei xL I ,YL(10),XR(10),YR(10) 


ANGLE OF INCLINATION IS BETWEEN 0 AND 89.9999 DEGREES. 
TAN(89.9999) - 572585000 
ТЕ ( ХМ .GE. O. .AND. Z2MOELTI 57295 GO ,ΠΠῚ 
DO 80. J =), Kal 


INY = YA J 
YL (MK) = = YA(J) + SQRT( LMIN**2 / ((1./YM **2)+1.)) 
INA = YL(MK) 


74 


IF ( INA .EQ. O ) GO TO 80 
M = XCORD( YM,INA,INY,INX ) 
YR(MK) = YA(J) - SQRT( RMIN**2 / ((1./YM **2)*41.)) 
INA = YR(MK) 
IF ( INA .EQ. O ) GO TO 80 
XR(MK) = XCORD( YM,INA,INY,INX ) 
MK = MK + 1 
80 CONTINUE 
x ANGLE OF INCLINATION IS BETWEEN 90.0001 AND 179.99 DEGREES. 
ELSE IF ( XM .LT. 0.) THEN 
DO 81Ј = 1, К, К-: 
INY = A 
INX = XA 
YL(MK) = D ο ο μι ο αι УМ *^*2)+1.)) 
ТМА = YL(M K) 
ТЕ (ТІМАВ.Е0. 04)%С60 TO 81 
H = XCORD( YM,INA,INY, ) 
YR(MK SN + SORT( ' AMINA / ((1./YM **2)+1.)) 
ТМА = IR (HK) 
IF ( INA .EQ. 0 ) GO TO 8l 
ME E YM,INA,INY,INX ) 
81 CONTINUE 
* ANGLE OF INCLINATION IS 90 DEGREES. 
ELSE IF ( ZM .GE. -572958.0) THEN 
Ρο 82 J= 1, K, K- 
INY = YA 
INX = ХА(Ј 
YL(MK) = a 
INA = L(MK 
IF ( INÀ т 0 ~) GO TO 82 
s = IN ІМІМ 
YR(MK) = YA(J) 
INA = YR(MK) 
IF ( INA .EQ. O ) GO το 82 
XR (MK) = INX + RMIN 
M MK + 1 
82 CONTINUE 
END IF 
RETURN 


END 
T T vedo eve eT v de Tee v vh e e e d vh e T cx 


* J. SUBROUTINE FOR COMPUTING THE RATE OF FLOW. 


T X X Zç Zc Zç e e T e dde e de de T6 T e e de de de de Te T T e 7k 7k Z Ze Ze Ze Ze 7 7k Zc Ze Ze Ze Zc 7 7 Zk Zc Ze ye Z Zc Z 7k Ze 7k dex 


SUBROUTINE FLOW(SP,WIDTH,ARCTP , CASE ,FR,MAXNO) 


INTEGER 

-  DIMENSION adt 
DIMENSION DW(3. Ἂς DD 
DATA DSTN  / 


АВСТР, СА5Е 
2),FR(S 
3, 2) 


5) ,DSTN(6) ,MAXNO(3,3) 
ο ου. 2. { 


T 
» DSTN IS DISTANCE BETWEEN EÁCH ELEMENTS. 


| I ARCTP 
* — ASSIGN DOCTRINAL WIDTH AND DEPTH OF BATTALION. 
W 


о 
= 
Q9 600 CO CO PO PO PO DB 2 КЗ 9S 


= 
н 
О 
τὶ 


` 


~ 


m 


~ 


= N) N) = i— N) N) У КР ЕЗ 
ЖЮ Ооо ооф оо 


Hl H il ill H H lH H H H H 


COMPUTE THE FLOW RATE. D BATTALION / HOUR ) 


IF ( ез ыс: | NI 


, 


~li tl 


р = DW(J,K) * DD(J, к) 
FR(J,K) = ( SP(I,J * WIDTH1 ) / D 
80 CONTINUE 
DO" 81 ο e 
IF(WIDTH .GE. 227 THEN 
WIDTH1 = 0.005 
END IF 
D= DAKE PES DD ү] AS. 
FR(J, K) = ( SP(I,J * WIDTH1 ) / D 
81 ONTINUE 
ELSE IF K κ. ЕО, 2” “ИНЕП 
50 δο. К 51/2 
IPK E 2 .AND. WIDTH1 .GE. 0.006) THEN 
W DTH1 = 0.006 
ELSE ч эше 2 .AND. WIDTH1 .LT. 0.006) THEN 
БАРЫН WIDTHl 
END 
D Dwg, ‚К) ж DD(J, M 
FR(J, К) = ( SP(I,. 2,5 * WIDTH1 ) / D 
83 CONTINUE 
END IF 
Χ COMPUTE THE МАХ όσον OF ELEMENT IN MULTIPLE COLUMN FORMATION. 
IF ( CASE .EQ. l ) THEN 
КМА = 
DO 85" ҚА - 1,2 
DO 85 KB= 1,2 
IF( KB .EQ. 1) THEN 
МАХМО (КА ,KB) = WIDTH / DSTN (KNA) 
ELSE IF(KB ΠΝ 
МАХМО KA RB) = (WIDTH1 1100" - DSTN(RNA) 
END IF 
KNA = KNA + 1 
85 CONTINUE 
ELSE IF ( CASE .EQ. 2 ) THEN 
KNA = 2 
DO 836 КВ = 1,2 
IF( KB .EQ. 1) THEN 
MAXNO(KA, 8) - D / DSTN(KNA) 
ELSE ТЕСКЕМЕ ОЯ 2 | 
МАХМО(КА, AE " WIDTH * 1000) / DSTN(KNA) 
END IF | | 
КМА = КМА + 1 
85 CONTINUE 
END IF 
RETURN 


END 
Xx h a ok na 


х Қ. SUBROUTINE FOR COMPUTING THE MIN TRAVERSAL TIME ALONG THE ARC.. 
kkkkkkkkkkkkkkkřkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkřř 
SUBROUTINE TIME(ARCTP,SP,UNITF,CASE,DIST,WIDTH,ATIME) 

INTEGER ARCTP , CASE , UNITF 

REAL KR 

DIMENSION Ы ооо, 

DATA KR / O. 
* KR IS COEFFICIENT OF ROUTE AVAILABILITY. 

1 = АКСТР 


76 


90 


25 


29 


BDIST 2 DIST / 1000; 
IF ( UNITF .EQ. 1 ) THEN 
IF( CASE EQ. οἱ » THEN 
DO 90 


E ADIST /( KR * SP(I,KU,1)) 
CONTINUE 
ELSE IF( CASE Е 1592,2 
35 2 ADIST РТ 3,1) 


END ТЕ 
ELSE IF m m EQ. 2 ) THEN 
τες C 6 1 ) THEN 
DO 25 Jus 
ATIME(J) = ADIST / SP(I,J,K) 
CONTINUE 
ELSE IF( m SEO; ) THEN 
185/57 = ADIST / SP(I,3,2) 
END IF 
END I 


F 
IF ARC WIDTH IS LESS THAN MIN FORMATION WIDTH, 
MAKE THE TIME BIG NUMBER. 
BIG 9.99 
FORWV 
FORWD = 2. 
IF ( UNITF .Е0. 1 ) THEN 
IF( CASE .ЕО. 1 ) THEN 
IF( WIDTH .LT. FORWV ) THEN 
Dô 97 KN = 1,2 
ATIME(KN) = BIG 
CONTINUE 


END 
ELSE S CASE .EQ. 2 ) THEN 
IF( WIDTH .LT. FORWD ) THEN 
ATIME(3) = BIG 
END IF 
END IF 
END IF 
RETURN 


END Р 
kkkkkkkkkkřkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


* L. SUBROUTINE FOR PRINTING THE RESULTS. 
eh T e e e e e he eh e e he e he e ehe e e e eh he e e he he e e he he e e hh e che e he e e he 


SUBROUTINE PRT 
INTEGER ΠΠ ΑκοτΕ ΕΕ CODE, PRINT ,UIMX ,UTMY, UNITE ,SPD 
REAL LMIN 


DIMENSION ALTD(100,100) ,NCODE(100,100 NX (10) ,NCODEX (10) 
DIMENSION XNODE (20) , YNODE(90) , CODE(0:7,0:7) 

DIMENSION DISTA(90),XA(90),YÀ(90) 

DIMENSION HT(90), DHT (20), ,SLOPE1 (90) 

DIMENSION DISTB(100,100 58100, 100), ,YB(100,100) 
DIMENSION HT2(100,100) , DHT2(100,100) , SLOPE2(100,100) 


DIMENSION NO 90), , SLPCH(90) 

DIMENSION  XL(10),YL(10),XR(10), 18 (10) 

DIMENSION SP 7,3,2),FR(5,5) ,MAXNO(3,3) 

DIMENSION ACTD(90) ,ATIME(3 

DIMENSION ROADW(7 

COMMON / INTS / NSTEP,IX1,I¥1,1X2,I¥2,JA,JB,K,NARC,NCOUNT, 
ж ΝΟΡΕΆ,ΝΟΡΕΒ, J1,J2,33,34, PRINT, UTMX, UTMY, 
ж CASE ,ARCTP,FR,AD, UNITF ,NSPCL 


COMMON / INTA / NCODE , NO, BLID, CODE, MAXNO 
COMMON / SREAL / XM, ANGLE, ,DIST, TSLPCH, ASLPCH,ALT1,ALT2,AVGSLP, 


* ІМІМ, ЕМІМ, AL, YL, АК, ҮҢ, M RY, WIDTH, THRES1, THRES2 
COMMON / AREAL / XA, ҮА HT, DHT, DISTA, DISTB, SLPCH, SLOPE1, SLOPE2, 

* AB, ,YB, ,HT2, DHT2, SP „АСТР, АТТМЕ, ROADW 

OPEN UNIT = 10, FILE = 'NETW1' 

OPEN UNIT = 20, FILE = 'NETW2' 
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* 


700 


701 
730 


τοι 
ӘЛ 


132 


8 
92 


734 


735 
92 

136 
740 


741 
742 
750 


του 
34 


752 


753 


754 


756 


OPEN UNIT 
OPEN 


29, FILE ' МЕТА! 


UNIT 26, FILE = ПЕТЕ; 


PRINT GENERAL INFO EE EACH ARC. 


жж ж 


+» 


IF ( PRINT EQ. 1 TH 

WRITE (6,700) COUNT , NODEA , NODEB 

FORMAT(! | ARC = пио 

/ 1, 4X. 1(!- Dl 71; 5X,'NODE ',12,6X,'NODE ',I2, 
X.' XM Lek ACT. DISTANCE', 

Е АХ,! 1Х1 ыл 5X,'IX2  192') 
WRITE (6, 701) IXl, IYl,IX2,IY2,XM,ANGLE,DIST,AD 
RMAT(, A ,1X,216,2X.216,1X,F842 PEM2 2X. F11.2) 2:08.02) 
BITE 730 ) 


FORMAT(' ',/,4X,71('-!'),//,20X, ' SOME POTTS ALONG  ARC!,/, 
4X, nem | I 'COORDINATE( X : Y )',5X,'ALTITUDE',/) 
DO 91 MA =1,K 


WRITE(6, 731) MA ,XA(MA),YA(MA) ,HT(MA) 

FORMAT ( ! Es 'POINT | ТЭТ О πα 2 врео 
EM 
WRITE(6,732) 
FORMAT N ‚Г.а, ПЕВ 
Τις, "POSITION ‚8Х 'DELTA H! „ох, ‘DISTANCE НБА SLOPES] 
‘ACTUAL DISTANCE! 
рО 92 MB = Ет 

WRITE(6, 733) МВ, MB+1 ,DHT (MB) , /DISTA(MB) , SLOPE1 (MB), ACTONI) 

FORMAT ( ' 5” IPOINT ' uo, το 212: ЕТАН 812. 1) 
CONTINUE 
WRITE(6, 734) 
ЕОКМАТ(' ',/, 7 
το, IPOSITION oM SLOPE CHANGE!) 
DO 33 MO S l RS 


PORC 235) MC,MC+1, SLPCH(MC) 


iS: 


RMAT( !,3Z, IPOINT ' ul37 213/2Х ΕΘ) 
СОМТІМ 
WRITE (6. , 738) TSLPCH,ASLPCH 
FORMAT ,48,' TOTAL SLOPE ος : ,10.4 


/,4Х, AVERAGE ‘SLOPE CHG. 3 I kO 4X. ЕКЕ 3) 
МЕТТЕ (6, 740) NODEA,NODEB 


FORMAT(' ',/,4X,'NÓDE' ,13, ': — ALTITUDE l ', 
4X, 'NODE', n. ADPITUDES2U 
WRITE(6,741) ALTI ALT? 
FORMAT(! ',20X,F5. i x, F5.1) 
WRITE( 745) AVGSLP 
FORMAT ,AX,'AVERAGE SLOPE ЕТ А Ақ пештен 
WRITE (6, 750 
FORMAT(' !',15X,' SOME POINTS OFF ARC ( FROM EACH POINT ', 
'ON THE are )',/,7X,'COORDINATE ( X Y )',8X,'ALTITUDE!, 
4X,' DELTA H',4X,'DISTANCE',5X, SLOPE',/) 
Ni = 0 
DO 94 JX =1,K 
N1 = NÍ+ 1 
DO 94 JY = 1,NO(N1) 
WRITE(6,751) JX,JY,XB(JX,JY), XS (o3 | HT2 (JK JY), 
DHT2(JX,JY),D TSTB (JH, JY) ,SLOPE2 J 
FORMAT(' ',3X,'POINT ',12,' 12. 287.1, ин)" АЕТ) 
CONTINUE 


WRITE (6,752) LMIN, RMIN WIDTH 
FORMAT ( ! |^ as TRE M 25X, MINIMUM DISTANCE OFF ARC', 
/,4X,'LEFT MIN DISTANCE ' :',F8.1, 
/,4X,'RIGHT MIN DISTANCE  :',F8.1 

/ 4X 'ARC WIDTH :',28.1) 

МЕТТЕ (6,753) 

FORMAT(' ',/,27X,'POINT 1',13Х,'РОІМТ 2', 
/,24X,'X LEFT Ү ГЕЕТ!,4Х,!Х RIGHT Y RIGHT!) 

WRITE(6, 754) (XL(K) ,YL(K), KR (IO ,YR(K), K-1,2) 
FORMAT(' ',/,4X,'INITIAL PO : ' 4F10.1 

"4X, ' END EGENT : ! 4F10.1 
WRITE (6,7 


RMAT(' ге), БАПАН υ η 
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ж 


760 


751 
762 


766 
25 


758 


26 


7620 


7660 
950 


7680 


МОЗ 


764 
65 
767 


Эб 
774 


5 
ж 


PRINT FLOW RATE AND TRAVERSAL TIME. 


a+ OE 


+» 


ELSE IF ( PRINT .E 2 ) THEN 
WRITE(20, 760) NCO T NODEA, NODEB 
FORMAT un ARC o la 
ах, 71('- E NL 5X, ' NODE ,12, ek. 'NODE les 
6X, ХМ ANGLE DISTANCE', 
А! ТУ I DIM ПА 
ЁТЕ (20, πο την το πο XM. ANGLE ή 
БМАТ(! 1,1Х,216,2Х,216,3Х,3Е11.2) 
RITE (20, , 762) LMIN, RHIN, ROADW(ARCTP), ЖЕРРИТАКСІР 


FO ,/ 4X, 71(' - ve 25Х, 'MINIMUM DISTANCE OFF АВЕ 
ЕК LEFT МІМ SISTANCE Бе, 

/,4X, ' RIGHT MIN DISTANCE mee, 

jee, WIDIH- OF-ROAD ITSELF :',F8.1, 

// ,4X,'ARC WIDTH е БӨРІ, 

Le 4X, АВС ΤΕΕ ΓΕΝ 

f, SOM UNTT ТУРЕ ,2Х, MISSION TYPE! , 3%, 


FLOW RATE',2X,' MAX £ ELEMENT!) 
IF ( CASE .EQ. 1) THEN 


Dos u-- I2 
WRITE(20,766) IJ,IK,FR(IJ,IK),MAXNO(IJ,IK) 
КОРШІ А, ІЗ 3/12, 103E7 . 2,9. L49 
CONTINUE 


ELSE IF ( E .EQ. 2) THEN 


PO ЗБЕ П 1 

НИК ВЕСЕО ЕК) МАХМОСТО, ІК) 
ПӘК 12-105 Е625,9Х,14) 
CONTINUE 


END IF 
WRITE(20,7620) 
EOE ',/,4X,'UNIT TYPE MIN TRAVERSAL TIME' ) 
РЕ САЗЕ с 1) ТНЕМ 
DO "950950 T2 
ΤΙ ΓΕ ATIME(UK) 
РОВМАТ Ὅς τ 9.59.2) 
CONTINUE 


ELSE IF ( CASE EO. 2) THEN 


WRITE(20, 7680) ІК,АТІМЕ(ІК) 
FORMAT(' ‚7Х,12,9Х,ЁЕ9.2) 


ІЕ 
WRITE (20, 763) 
μι 7 27 τοπ) 13 ШРОТМТ 2', 
ы 24X,'X LEFT Y LEFT! 4X, 'X RIGHT У RIGHT!) 
НАІТЕ (20, 764) ЖЩ), YL(1).XR(1),? ?Е(1 
FORMAT /,4X,' INITIAL POINT : ',4F10.1) 
RITE (20, 783) uc. 212). ER), 'vR(2) 


,4210.1 
WRITE(20, 767 ЈЕ 
FORMAT(' ',/,4Х,71('*')) 


PRINT COORDINATES OF BOUNDARY LINE. 


SESH ae ot PRINT p 4 ) THEN 


DO 97 IN s 1, 
XL(IN) = Utitx +(XL(IN) /1000. 
YL(IN) = UTMY +(YL(IN) /1000. 
XR(IN) = UTMX +(XR(IN) /1000. 
YR(IN) = UTMY +(YR(IN) /1000. 
CONTINUE 


WRITE(10,774) NODEA, XL(1),YL(1),XR(1) УК(1) 
ВМАТ(! !,/,4%, 

МЕТТЕ (10, 75) τι x. КИ 0025 YR(2) 

ҒОБМАТ(! 1,3%, !МОПрЕ!,12 ! AF10.2) 


PRINT NETWORK DATA FOR CARTESIAN SPACE NETWORK PROGRAM. 


BLE IF ( PRINT .EQ. 5 ) THEN 
IF( NODEA .GT. NODEB ) THEN 
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874 


NODE = NODEB 


NODE2 = NODEA 

ELSE ТЕС ^ NODEB .GT. NODEA ) THEN 
МОРЕ1 = МОРЕА 
NODE2 = NODEB 


END IF 
DST = DIS? / 1800: 
WD = WIDTH / 1000. 


IF( CASE .50. 1 ) THEN 
Кб = 


KT | 
ELSE IF( == EOL 20) TEN 
KS 


KT 
END IF 
NTIME = ATIME(KS) * 100 
NFR = FR(KS,KT) * 100 
TARCT =) Shere 


ITF 

SPD = SP( IARCH κ, IUNITEM 
Το CASE EO NEED THEN 

WRITE (28, 874) NCOUNT ,NODE1 ,NODE2 ,NTIME,NFR,DST,WD,SPD 
ELSE ТЕ(БСАЗБ .EQ. 2 THEN 

WRITE ey e: NCOUNT,XODEl1,NODE2,NTIME,NFR,DST,WD,SPD 


FORMAT(' '. .,1X,13,215,216,3X,F4.1,3X,F6.3,2X, I4) 
RETURN 
END 
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COMPUTER 


APPENDIX D 
PROGRAM FOR CARTESIAN SPACE NETWORK 


This appendix contains the computer program used to support the algorithms in 


Chapter IV. 
PROGRAM NET 
ж JULY 13, 1987 (13:00) 
A NETWORK MAIN PROGRAM 
* SHORTEST PATH ALGORITHM 
* Xx X Z Z Z Zç Zc Yc Zc Zc Z κκ VARIABLE DECLARATION ας Ἂς ας Z Z e ve e de Ж 
INTEGER TAIL,HEAD , TIME , FLOW 
INTEGER CASE,NTYPE,PARAM,SP,UNITF,ROWSP 
CHARACTER *1 RESM 
DIMENSION 


DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENS ION 


COMMON / INTS 


COMMON / 


COMMON / 


E. COMMON / 


FLOW(100) ,DIST(100) ,WIDTH(100) 
LINK(100 , 100) LARC (100) 

NO 100) SP(100) 

INODE(1$0} 3 JNODE(150) , ARCOST(150) 


І5РАТН(5 
CASE ,NTYPE,N,M,NARC ,ISTART ,LAST, PARAM, 
NUMP , UNITF ,NVEH,NCOL,NROW,ROWSP ,MINE, 
NWD , NDP , NFORWD 

INTA / TAIL,HEAD,TIME,FLOW,LINK,LARC,NO,SP, 
INODE , JNODE , ISPATH 

REALS/ BIG,TPATH,XLEN ,PATHL ,FORWD,FORDP,PRT 

REALA/ WIDTH,DIST,RESM 


ΠΗ ЁТТМЕ( 150.) 


к (1). SET THE INITIAL OPTIONS. 





CALL OPTION (CASE,UNITF,NCOL,ROWSP,PARAM,FORWD,NFORWD, 


RESM,MINE) 
SET THE START NODE 


ISTART = 1 
сы SET THE LAST NODE 
LAST = 34 
es TOTAL NUMBER OF ARC IN THE SECTOR 
NARC = 72 
s TOTAL NUMBER OF NODE IN THE SECTOR 


ы SET THE INITIAL CONDITIONS. 


CALL INIT 
ж 


(ARCOST,NUMP,ISPATH,XLEN,ADJUST,PATHL,SPEED, 
LINK,TPATH,NNODE,BIG) 





Б (2). READ THE DATA FROM THE DATA FILE.(DATA 1, OR DATA 2 ) 





x DATA 


1 ; ARCS AND CHARACTERISTICS FOR VEHICLE UNIT 


^ DATA 2 ; ARCS AND CHARACTERISTICS FOR DISMOUNTED TROOPS 
_ DETERMINE VARIABLE FOR READING INPUT DATA.( VEHICLE OR DISMOUNTED ) 
IF( PARAM .Е0. 1 THEN 


IF( CASE .EQ. 1 ) THEN 
KUNIT = 

ELSE IF( CASE .EQ. 2 ) THEN 
KUNIT = 2 

END IF 


ELSE ae ld e „EO. 2) THEN 


END 
READ(KUNIT, loo, END= =999) (wo (I ^ ,TAIL(I),HEAD(I 
OW(I),DIST(I),WIDTH(I),SP(I) 
100 FORMAT (2X,I13, 215. 216. 3X,F4.1,3X,F6.3,2X, 14) 


* SPECIFY THE WIDTH AND DEPTH OF UNIT FORMATION. 
2358 CALL  FORMIN (NARC,CASE,UNITF,NCOL,ROWSP,PARAM,WIDTH,TIME, 
* NROW,DIST,FLOW,FORWD,FORDP,NVEH,NWD,NDP) 


), TIME(I) 
, I21,NARC) 


* 


а (3). DETERMINE THE APPROPRIATE VALUE FOR PARAMETER. 


* DETERMINE THE VALUE OF PARAMETER WHEN THE MINEFIELD IS ON ARC. 
IF Il .EQ. 'Y') THEN 
F (PABAM .EQ. 1) THEN 
TIME(MINE) - TIME(MINE)* 4 


END 
END IF 
* DETERMINE THE APPROPRIATE VALUE FOR PARAMETER 
N = NNODE 
M = NARC 


DO 10 KA = 1,NARC 
Ὃ вен ш ος 
JNODE(KA) = HEAD(KA 
THEN 


IF (PARAM EQ. 1) 
) = TIME(KA) 
ELSE IF CPARAM -EQ. 2) THEN 


OST(KA) = DIST(KA) 
END 


10 CONTINUE 


s (4). DETERMINE THE MIN (PARAMETER) PATH. 


„CALL SHORTP (N,M,INODE,JNODE,ARCOST,ISTART,LAST,BIG,NUMP, 
ISPATH, XLEN, NP ) 


* DETERMINE THE ARC NUMBER ALONG THE PATH. 
DO 60 М3 NARC 
LINK(TAIL (M3). HEAD(M3)) = NO(M3) 
60 CONTINUE 
PATHL = 0. 
DO 70 MA = 1, NUMP-1 ы 
LARC(MA) = LINK(ISPATH(MA), ISPATH(MA*1)) 
PATHL = PATHL * DIST( LARC(MA) ) 
70 CONTINUE 


* WHEN 'PRT' IS 1.0, THERE 15 159 FEASIBLE ROUTE TO GO. 
ТЕ (РА RAM EQ. ET 
РЕТ = 0. 
DO 75 M MB = 
IF (TIME (LARC(MB) ) .GE. 999) THEN 


END IF. 
75 CONTINUE 
END IF 
х DETERMINE THE s SEED = O. TIME FROM START NODE TO LAST NODE. 
E 


DO 80 M5 = NUMP 
SPEED = SPEED + (DIST (MS) * SP(M5) / PATHL ) 


80 CONTINUE 
ADJUST FORDP / SPEED 
XLEN + ( ADJUST * 100. ) 


H H 


TPATH 


82 


т (5). РКІМТ THE RESULT ОЕ CALCULATION. 


CALL PRNT 
SLOP 


END 
kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


х А. SUBROUTINE FOR DETERMINING THE INITIAL OPTIONS. 


ΧΧΧ — Zc Zc Z< Zç Z< ç % > ç Z Z< >< e Te ee Too de Zc Zk 7 Zç Zç > Z Zç Zc Zç Zç k Zç Sk Ze 9 7 ç Zc ded %—< Z *—< =< 


SUBROUTINE OPTION (CASE,UNITF,NCOL,ROWSP, PARAM, FORWD ,NFORWD, 
ж RESM, MINE) 

INTEGER CASE, PARAM, UNITF, ROWSP 

CHARACTER  *1 RESM 


WRITE(6, 603) 


503 FORMAT(! !,///) 

ШЕ PRINT* 
PRINT*, | WHICH TYPE OF PARAMETER DO YOU DESIRE?! 

I 

EA ) X fece Reef d e defe e sehe ee eee 
PRINT*, | ж τι 
PRINT*, ' * Га MIN TIME PATH ж! 
PRINT*, | ж ж |! 
PRINT | ^ 2. MIN DISTANCE PATH x! 
PRINT*, ' < жа 
РЕТМТХ, | d Yee sese dee ee ehe seo ie Je ese eee dci e eee 


READ(5,*) PARAM 
PRINT*, '** NOTE : FOR YOUR REFERENCE, CURRENT ANSWER IS', 


* PARAM 
IF ( PARAM .EQ. 1 ) THEN 
GO TO 16 
ELSE IF ( PARAM .EQ. 2 ) THEN 
со TO 20 . 
ELSE IF ( PARAM .NE. 1 .OR. PARAM .NE. 2 ) THEN 
PRINT*, ' 
PRINT^*, '*** ERROR : ENTER THE NUMBER 1 OR 2 ***! 
PRINT*, | 
GO TO 15 
END. ТЕ 
16 WRITE(6,601) 
601 FORMAT η 
PRI! МІХ, | WHICH TYPE OF UNIT DO YOU DESIRE? ( 1 OR 2 )' 
PRINTS, RM νν -е 
PRINT*, | * жі 
PRINT л КООУЕНТСГЕ UNIT τ. 
PRINT ' * x i 
PRINT^ ' * 2. DISMOUNTED TROOPS zd 
PRINT*, ' * τ 
PRINT*, ' ЖАКА 
READ(5,*) CASE 
x PRINT*, '** NOTE : FOR YOUR REFERENCE, CURRENT ANSWER IS', 
CASE 
IE CASES .S0w-L .ORGCASE .EO& 2 ) THEN 


1 
ELSE DELI CASE КИЕ КОЛОН. CASE .NE. 29) НЕМ 


PRINT” 
PRINT^, ο ERROR : ENTER THE NUMBER 1 OR 2 ***' 
PRINT*, | | 
GO TO 16 
Ер. m 
17 WRITE(6,602) 
602 FORMAT (' S 
PRINT^, ' WHICH TYPE OF FORMATION DO YOU DESIRE? ' 
PRINT ! ee RA AO RIEN 
PRINT*| I x κ |! 
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πα Ξ 1. MULTIPLE COLUMN FORMATION de 
PRINT*, ж 2. SINGLE COLUMN FORMATION wl 
PRINTS, * κι 
PRINT*, ! cce hehe heh e eee hehehe hh he he he he heh e hee eee 


READ(S5, ,*) UNITF 
PRINT*, '** NOTE : FOR YOUR REFERENCE, CURRENT ANSWER IS', 


UNITF 
IF ( UNITF .EQ. 1 .OR. UNITF .EQ. 2 ) THEN 
GO TO 18 


ELSE IF ( UNITE «ΝΕ. 1 .OR. UNITF .NE. 2 ) THEN 
PRINT*, ' 
PRINT*, ххх ERROR : ENTER THE NUMBER 1 ОН 2 ХХХ! 
PRINT*' : 
GO A 17 
END 
18 IF (CASE Pd T TH 
UNITE "EO. 1) THEN 
WRITE(E ,60 24) 
604 FORMAT ( ' РАДИ ο ή 
PRINT ENTER THE NUMBER OF COLUMN.(5, 10,.20)' 
D. 
PRINT*, i * * * * * * ж! 
PRINT* l i 
PRINTZ, i ж * * * * x ж! 
PRINTS, 
PRINT* i * κ * * * * ж 
READ(5, 2» NCOL 
PRINT* 
WRITE (6, 6 05) 
605 RMAT (' v 
PRINTS 17 | | ЕМТЕН THE DISTANCE BETWEEN EACH ROW. (75550 ‚25)' 
PRIN 
Ати i / * * * * * * ж! 
РЕТМТА, ! || | 
PRINT*, i / * * * * * * ж! 
ΡΗΙΝΤΑ, ' : 
PRINT*, i * * * * * * ж! 
РВФМТ ы! ' i 


READ(5, ,*) ROWSP 
ELSE IF (UNITF .EQ. 2) THEN 


WRITE(6,60 96) 
606 FORMAT (' ΜΙ: 
PRINT*, | ENTER THE DISTANCE BETWEEN EACH ROW.(75,50,25)' 
аы. i / * * * * * * жі 
PRINT*, | || i 
PRINT*, ! / ж * * * * * жі 
PRINT, | | | 
PRINT*, i * * Ж * * * жі 
РЕТКЕ i)! 
ΕΕ ж) ROWSP 
ND IF 
ЕБЗЕ ІІБ (CASE с 2) ТНЕМ 
UNITÉ .EQ. 1) THEN 
WRITE (6,2 14) 
614 FORMAT (' LL i 
| PRINT*, ' ENTER THE NUMBER OF COLUMN.(10, 20, 50)! 
PRINT*| | 
PRINT*, ! х * * * * * *! 
PRINIC ! 
PRINT*, i * * * * * * жі 
ES | 
PRIN ! * * * * * * Ж! 
READ( 5, ‚жу м NCOL 
PRINT <a 
WRITE(6,615) 
615 FORMAG ТТТ O) 
ἘΠ i | ENTER THE DISTANCE BETWEEN EACH ROW.(5,10,15)' 
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ды ο πο ο ~ хл әх! 

КІ 

PRINT^, ' | ko ok х xk k ж x! 

PRINT*, | 
4 


PRINT, * * * * * * ЖІ 


PRINT ж, 
EE ‚х ROWSP 
END IF 


END IF 
20 WRITE (6,607) 
607 FORMAT (' E) 
РЕТМТХ, IS THERE ANY MINEFIELD ALONG THE ARC? (Y OR N)' 
READ(S,' (aL) ) RESM 
IF (RESM .EQ. 'Y') THEN 
PRINT*, ! ENTER THE ARC NUMBER OF MINEFIELD ARC' 
READ(S,*) MINE 
ELSE IF (RESM .EQ. 'N') THEN 
MINE = 0 
END IF 
RETURN 


END 
kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


* B. SUBROUTINE FOR DETERMINING THE INITIAL CONDITIONS. 


kkkkkkkkkkřkřkkkkŘkřkkkkŘkkkkkkkkkkkkkkkkkkkŘkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


„SUBROUTINE INIT (ARCOST ,NUMP , ISPATH,XLEN ,ADJUST, PATHL,SPEED, 
LINK, TPATH, NNODE Резо) 


DIMENSION ISPATH(50) ,ARCOST(100) 
DIMENSION LINK(100,1 100 0) 


ВІС 1.0Е10 
NUMP = 
XLEN = 
ADJUST 
PATHL 
SPEED 
TPATH 
DO 10 


10 CONTINUE 


PE II li Ji ооп 
" I| OO O O 


20 CONTINUE 
DO 207 η]. 
DO О Ја 
LINK(M1,M 
η CONTINUE 
RETURN 


END 
kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkěkk 


ІШІ 


х С, SUBROUTINE FOR DETERMINING THE FORMATION WIDTH AND DEPTH. 


„SUBROUTINE FORMIN (NARC,CASE,UNITF,NCOL,ROWSP,PARAM,WIDTH,TIME, 
NROW,DIST,FLOW, FORWD, FORDP, ,NVEH, NWD ,NDP) 


INTEGER TIME,FLOW 
INTEGER CASE , PARAM, UNITF ,ROWSP 


DIMENSION TIME(100),FLOW(100) ,DIST(100) ,WIDTH(100) 
8161 = 999. 


х 


T (A). DETERMINE THE FORMATION WIDTH ОҒ А UNIT OF MOVEMENT. 
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ТЕ( CASE .EQ. 1 ) THEN 
NVEH = 50 
IF ( UNITF EQ. 1 ) THEN 
NROW - NVEH / ΝΕΟΙ. 
COLSP = 25. gl 
FORWD = A (NCOL- 1) * COLSP) / 1000 
ELSE τες UNIT ТЕ .EQ. 2 ) THEN 


NCOL - ΝΤ / NROW 


С) 
О 
t 
оз 
“Ὁ 
II ll 


Es FORWD (NCOL * COLSP) / 1000 
FORDP = ((NROW-1) * ROWSP) / 1000 
ELSE IF( CASE .EQ. 2 ) THEN 
ж WHEN UNIT IS MULTIPLE COLUMN FORMATION, THE FORMATION WIDTH 
ж WAS ENTERED. 


IF ( UNITF EQ. οἱ ) ТНЕМ 
ММЕМ - 
NROW = AMEN / NCOL 


COLSP - 
FORWD NCOL * ROWSE, 
FORDP NROW .* ROWSP 
ELSE IF (UNITE .Е0- 2 ) τα 
FORWD = 0.006 i 
FORDP = 2.5 
ЕВЕ 
NWD = FORWD * 1000 
NDP = FORDP * 1000 


/ 1000. 
/ 1000. 
N 


END IF 


А (В). ASSIGN BIG NUMBER WHEN ARC WIDTH IS LESS THAN FORMATION WIDTH. 


DO 50 І 1, 
IF (PARAM . „EQ. D THEN 
IF (WID a . LT. FORM ТНЕМ 
МЕ(І) = BIG1 
END 
END IF 
50 CONTINUE 
RETURN 
END 


ХХ e ЗОО ЈОЈ О О О О О О do dx 


* D. SUBROUTINE FOR DETERMINING THE SHORTEST PARAMETER PATH. 
ΑΛΑΛΑΛΑΛΑΛΑΛΑΛΑΛΛΑΛΑΛΑΛΑΛΑΛΑΛΑΛΑΛΑΛΑΛΑΛΑΛΑΛΑΛΑΛΑΛΑΛΑΛΑΛΑΛΑΚΑΛΛΑΛΑΛΑΛΑΛΑΛΑΧΑ 


,SUBROUTINE SHORTP (N, M,INODE,JNODE,ARCOST,ISTART,LAST,BIG,NUMP, 
ISPÁTH,XLÉN,NP ) 

х MINIMUM PATH 

x FIND A SHORTEST PATH BETWEEN TWO GIVEN NODES 


INTEGER INODE(M) ,JNODE(M), ISPATH(N) 
REAL WK4 , ARCOST(M) 
LOGICAL ORE ,IFIN 
DIMENSION IWORK2(34) , INORK3(34), WK4 (34) 
ж 34 : NUMBER OF NODES IN NÉTWO 
IWORK2 


DO 10 Is1,N 
ЧКА (1) 
I) 
10 CONTINUE 
WK4 (ISTART) - 0. 


- ΙΙ И 


‚ТЕШЕ. 
0 


EU 
= ISTART 
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IWORK1 (ISTART) = .FALSE. 
МР = 0 


XLEN = Q. 
х 





= (A) FOR EACH FORWARD ARC ORIGINATING AT NODE I CALCULATE THE LENGTH 
OF THE PATH TO NODE I. 


Í 





20 Ic = 0 
DO 30 K Ded 
TE (коо (К) E 
= IC + 
po di 
ΤΕ 


І) ТНЕМ 


li du cuo 


ISPATH 


END IF 
ТЕ ANODE (K) 
= IC + 
monka IC) 
ISPATH(IC 
END IF 
30 CONTINUE 
ЈЕ CIC беле 49 T THEN 
DO 0 РС 
товка (Е, 
ISPATH(L 


(INORK1(J)) THEN 
z нка(т) + ARCOST dy À 
F (D .LT. WK4 THEN 
D 


WKA(J) s . 
CIE zu 
END IF 
END IF 
40 CONTINUE 
END IF 


K 
JNODE (K) 
I) THEN 


tO 


H Hit 


K 
INODE(K) 


K 
J 
ТЕ 
р 


х (В). FIND THE MINIMUM POTENTIAL. 


D = BIG 
IENT = 0 
IFIN = .FALSE. 
DO 50 I = 1 
Та (тйовкі(т)) ТНЕМ 


IFIN = UE. 
ΤΕ (WKA (1) JD D) THEN 
D = WK4(I) 
IENT = I 
END IF 
END IF 
50 CONTINUE 


х (су. INCLUDE THE NODE IN THE CURRENT PATH. 
x 


IF (D .LT. BIG) THEN 
IWORK1(IENT)=. FALSE. 
IF (IENT .NE. LAST) THEN 
I = IENT 


IF (IFIN) THEN 
NP = 1 


RETURN 
END IF 


END IF 


IJ - LAST 
UMP = 1 


N 
ISPATH(1) 


60 K = 


E 


T 
L 
J 
70 

Қ = 


00 


ISPATH 


70 J 
XLEN т 
RETURN 
END 

JC ke oe de 


с 


Е. SUBROUTINE 


W Z Z Z 7 Ze Zk Z< Ze Zc Zç 2k κκ 


ж 


SUBROUTINE 


INTEGER 
INTEGER 
CHARACTER 
DIHENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
COMMON / 


жж 


COMMON / 


COMMON / 
COMMON / 
OPEN cee 


OPEN 
* (A). PRINT THE 
Жж 


TER 
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ELS 


жж 


UNIT = 


= LAST 
IWORK2(IJ) 
IF (INODE(K) E 
IJ= JNODE(K) 
ELSE 
= INODE(K) 
END IF 


NUMP = NUMP + 1 
ISPATH(NUMP) IJ 
(τα МЕ. І5ТАҚШ 00 TO 60 
NUMP / 2 | 
NUMP 


IJ) THEN 


Т 
ISPA 


( 


1 5 
TH(I 


T 


) 
ISPATH(J) 
J K 


i 
WK4 (LAST) 


kkkkkkkkřřkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


FOR PRINTING THE RESULTS OF CALCULATION. 
KKAKKAKAAKAKKAAAKKAAKKAKKKKAAKKEKKAAKKRAERAKAAKKKKAAKKARK 


PRNT 


TAIL, HEAD TINET T LOW 
CASE , NTYPE , PARAM, SP , UNITF , ROWSP 
^1 ВЕЗИ 


ШАП гоа аа, ,ТІМЕ(150) 
FLOW(100),DIST(100),WIDTH(100) 
LINK(100,100),LARC(100) 

М0(100),5Р(100 

ENODE (150), , JNODE (150) , ARCOST (150) 

ISPATH(50 

INTS / CASE,NTYPE,N,M,NARC,ISTART,LAST,PARAM, 
NUMP,UNITF,NVEH,NCOL,NROW,ROWSP,MINE, 
NWD,NDP,NFORWD 
TAIL,HEÁD,TIME,FLOW,LINK,LARC,NO,SP, 
INODE, JNODE,ISPATH 

REALS/ BIG,TPATH,XLEN,PATHL,FORWD,FORDP,PRT 
REALA/ WIDTH, DIST, RESM 


= 10, FILE TIME! 
20 FICE ШАА 


ІМТА / 


RESULTS ОГЕШ ΤΙΣ ΕΠΤΗ 


У TIME PATH',/,5X,15('=-')) 
NROW,NCOL, ROWSP 

NROW , NCOL , ROWSP 
X,'VEHICLE UNIT ',/,5X,'FORMATION :',I3, 
,13,' COLUMNS',/, 17X,'ROW SPACE', 

I3,' METERS !) 

IF (CASE .EQ. 2) THÉN 
NWD = FORWD * 1000 
NDP = FORDP * 1000 
IF (UNITF .EQ. 1) THEN 
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tj 
tL 


E 


713 


714 


121 


722 


WRITE ( 713) NROW , NCOL , NWD , NDP 
WRITE(10,713) NROW,NCOL,NWD,NDP | 
КОВМАТ(' ',4Х, 'DISMOUNTED TROOPS | 1, 5%, 
' FORMATION шыка! Εις ντο... COLUMNS ! Пе /, 
16X,' WIDTH',I5, E 15,' МЕТЕР5' | 
ELSE T (UNITF .EQ. 2) THEN 
WRITE ( 6, 21) NWD,NDP 
WRITE(10,714) NWD,NDP 
FORMAT(' ',4X,'DISMOUNTED TROOPS! ,/,5X, 
' FORMATION :',' WIDTH',I5,', DEPTH'.I6,' METERS') 
END IF 
END IF 
IF (RESM a 'Y!) THEN 


WRITE(10,721 
FORMAT (/, 55,11 MINE FIELD ( ARC NUMBER :',I3,')' 


ELSE IF (RESM o πρ ) THEN 
WRITE( 6,722 
ЕТТЕ 10. 722 
FORMAT(/,5X,'NO MINE FIELD ARC!,/,5X,35('-')) 


WRITE 71) MEE 


END IF 


* WHEN PRT EQUALS 0.0, THERE Y TH μα КООТЕ ТО σο. 
Е 


131 


732 


733 


734 


735 


736 


750 
760 


* 


IF (PRT шт 0.0 
ІҮ = 100 


ISUMT = XLEN / 100 
NXLEN = XLEN 
MSUMT = ( MOD( NXLEN, IY ) * 60 ) / 100 
ee HER ISUMT, MSUNT 
WRITE(10,731 ISUMT, MSUMT 
FORMATO UT / SK, 'SUM OF TRAVERSAL TIME : ! 
ПАРОХ НОК SSD S 2X. "MINUTE | » 


ITHR TPATH / 100 
NTPATH ТРАТН 
MINUT ( MOD( NTPATH, IY ) * 60 ) / 100 


WRITE 5,732) ITHR, MINUT 
ITHR. MINUT 


πο ',/, 5% 'TOTAL TRAVERSAL TIME 


Iu, HOUR',lI15,2Z, MINUTE) 


WRITE ( 6,733) PATHL 
WRITE(10,733) РАТНЕ 
FORMAT(' !',/,5X,'TOTAL PATH LENGTH ΤΙ. 
22k, KM! ) 
НЕДЕ, 6, dm 
WRITE(10,734 
FORMAT(! ',/,5X,'NODE NUMBER ALONG MIN TIME PATH ') 
WRITE | ЕНІН mu I =1, NUMP 
WRITE(10,735)(I,ISPATH(I), I зі, NUMP 
FORMAT (10X,13, 3X,13) 
WRITE ЖЕЗ 
WRITE(10,736 
ОВМАТ(! '!,/,5Х,! ARC NUMBER ALONG MIN TIME PATH ') 
WRITE | 6,730} (1 ,LARC(T), I -1, NUMP-1 
WRITE(10,750)(I,LARC(I), I =1, МОМР-1 
ОВМАТ(10Х,13,3Х,13 
WRITE(6,760) XLEN 
FORMAT ( | Шығыны. NUMBER ALONG MIN TIME РАТН :', 


* WHEN PRT EQUALS 1.0, THERE IS NO FEASIBLE ROUTE TO GO. 
ЕШӘЕ ТЕ (PREP EQ. 1.0) THEN 


741 
* 


š (В). 


WRITE(10,741 
FORMAT ( ,'THERE IS NO FEASIBLE ROUTE TO GO.') 
END IF 


2 
E 6,741 


PRINT THE RESULTS FOR HINIMUM DISTANCE PATH. 
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ELSE IF (PARAM .EQ. 2) THEN 
WRITE( 6,811 
WRITE(20,811 
811 FORMAT(5X,35('-'),/,5X,'MIN DISTANCE PATH',/,5X,19('-')) 
IF (RESM .ЕО. '4') ТНЕМ 
WRITE ( 6,821) MINE 
WRITE(20,821) MINE 
821 FORMAT (/ 5X, '1 ЊЕ Sp ( ARC NUMBER :',I3,')! 


/, 5 
ELSE IF (RESM .EQ. 'N' ) THEN 
НЕШЕ, Б 822 
WRITE(20 , 822 
822 FORMAT(/,5X,'NO MINE FIELD ARC',/,5X,35('=-')) 


END IF 
WRITE ( 6,831 
WRITE (20,831 

831 FORMAT(' !',/,5X,'NODE NUMBER ALONG MIN DISTANCE PATH ') 
WRITE( 6, 88: ISPATHIT], I =1, NUMP 
WRITE(20,832)(I,ISPATH(I), I =1, NUMP 

832 FORMAT (10X,13,3X,13) 
WRITE | 6,833 
WRITE (20,833 

833 FORMAT(' ',/ 5X,' ARC NUMBER ALONG MIN DISTANCE PATH ') 
WRITE | Ена DARET. I z1, NUMP-1 
WRITE (20,834)(I,LARC(I), I =1, NUMP-1 

834 FORMAT (10X,13,3X, 13 
WRITE | 6,360) XLEN 
WRITE (20,860) XLEN 

860 FORMAT(' ',/, SK, TOTAL DISTANCE ALONG MIN DIST PATH :', 


END IF 


* PRINT THE INPUT DATA FOR CHECKING. 
WRITE(S, 921) 
МЕТТЕ (6,920 (NO(Z) TAIL(I),HEAD(I),TIME(I),FLOW(I), 
ST(I),WIDTH(I),SP(I),I=1,NARC) 
IF (PARAM EQ. 1) THEN 


0 
ELSE IF (PARAM .Е0. 2) ТНЕМ 
MU = 20 


END IF 
ο ρα 921) 
WRITE(MU,920) (NO(I),TAIL(I), МЕР иар 
DIST(Ij, WIDTH(I),SP(I),IZ1,NARC) 
221 ЕОБНАТДЫ. 1270 ΥΣ МО TAIL HEAD TIME FLOW DISE 
< ' WIDTH SPEED! 
920 FORMAT (2X, 13,215 ,216pon be.) ok ho = се 1 
RETURN 


END 
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APPENDIX E 
COMPUTER EXEC PROGRAM 


This appendix contains two computer exec programs used to support the 
algorithms in Chapter III and IV. 
SINGLE ARC EXEC PROGRAM 

This exec program is used for the computer program for single arc attributes. 


&TRACE OFF | 

&TYPE Please provide the FILENAME for your VS FORTRAN program. 
&READ VAR &FN | 

&IYPE Do you need to compile your program ? (Y) 

&READ VAR &R COMPILE. . 

A SR COMPILE NE Y &GOTO -RUN 

БЛИШОӘКТУЗ ХЕМ 

ERC 20 0 &SKIP 9 

&TYPE Your program did not compile; check for errors. 
&IYPE Do you wish to XEDIT the program zile? (Y) 

&READ VAR &RESPI 

ОСИ СЕБРІ ПЕ У СЕХЕТ 1 

&COMMAND XEDIT &FN FORTRAN À | 

“ТҮРЕ Do you wish to run the program again? (Y) 

&READ VAR &RESP2 

EDSEXRESPA2CEO- Y &GOTO -H 

&EXIT 1 | | 
-RUN «ТҮРЕ Do you wish your INPUT to be from the terminal? (Y) 
&READ VAR &IN 

ΙΝ NE Y «GOTO -RUNZ . 

SUPUT FILE FILEDEF Ql DISK NB71 DATA Al 

ШЕШ” CILE FILEDEF 02 DISK NBNODI DATA Al 

ӘШШЕШТ FILE FILEDEF 03 DISK NBARC1 DATA Al 

-INPUT FILE FILEDEF 04 DISK SPEED РАТА Al 

FILEDEF 05 TERMINAL | ‚ 
-RUN2 &TYPE Do you wish your OUTPUT to go to the terminal? (Y) 
&READ VAR &QUT 

БЕН OUT NE Y &GOTO -OUTPUT ΕΙΤΕ 

&GOTO -LOAD 

“OUTPUT_FILE FILEDEF 06 DISK ХЕМ OUTPUT A (LRECL 133 
-LOAD LOAD &FN (START 

&iF XRC EQ 0 &SKIP 9. 

Ens Your program did not run correctly; check for errors. 
&TYPE Do you wish to XEDIT the program file? (Y 

&READ VAR &RESP3 

RERSUERESPS NE Y &BEXIT 2 

&COMMAND XEDIT &FN FORTRAN A 

&TYPE Do vou wish to run the program again? (Y) 

&READ VAR &RESP4 

&IF &RESP4 EQ Y &GOTO -H 

GBXIT 2 

CIPR OUT 30-Y &GOTO -REDO | 

СЕЕ τους output is in the file &FN OUTPUT A 

&TYPE Do you wish to BROWSE your output? (Y) 

&READ VAR &RESP 

&IF &RESP EQ Y &COMMAND BROWSE &FN OUTPUT A 

&TYPE Print your output file? (Ү) 

&READ VAR &RESP7 

ao EQ Y &COMMAND PRINT &FN OUTPUT A 


&TYPE Do you wish to XEDIT the program file? (Y/N) 
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&READ VAR &RESP5 

&IF &RESP5 EQ Y XEDIT &FN FORTRAN A | 

&TYPE Do you wish to run the program again? (Y) 
&READ VAR &RESP6 

&IF &RESP6 EQ Y AND &RESP5 EQ Y &GOTO -H 

p E EQ Y &GOTO -RUN 


NETWORK EXEC PROGRAM 
This exec program is used for computer program for Cartesian.space network. 


&TRACE OFF | 

“ТҮРЕ Please provide the FILENAME for your VS FORTRAN. program. 
&READ VAR &FN | 

“ТҮРЕ Do you need to compile your program ? (Y) 

&READ VAR &R COMPILE 

&IF &R COMPILE NE Y &GOTO -RUN 

-Н ЕОКТУ5 ХЕМ 

&IF SRC EQ 0 &SKIP 9 | 

“ТҮРЕ Your program did not compile; check for errors. 
&TYPE Do you wish to XEDIT the program file? 

&READ VAR &RESPI 

СТЕ SRESP! NES vl 

&COMMAND XEDIT ΔΕΝ FORTRAN A | 

“ТҮРЕ Do vou wish to run the program again? (Y) 

&READ VAR &RESPZ2 

&IF &RESP2 EQ Y &GOTO -H 

&EXIT 1 | | | 
-RUN &TYPE Do you wish your INPUT to be-from the terminal? (Y) 
&READ VAR &IN 

&IF &IN NE Y &GOTO -RUNZ 

-INPUT FILE FILEDEF O1 DISK МЕТА DATA А1 

TENPUT RILE ЕБІБЕРЕЕ 02 DISSEIIEIB DATA Al 

FILEDEF 05 TERMINAL | | 
-RUN2 “ТҮРЕ Do you wish your OUTPUT to go to the terminal? (T) 
&READ VAR &OUT 

&IF &OUT NE Y &GOTO -OUTPUT "ILE 

&GOTO -LOAD : 

-OUTPUT. FILE FILEDEF 06 DISK &FN OUTPUT A (LRECL 133 
-LOAD LOAD &FN (START i 

«ТЕ @RC БООР БКО 

&TYPE Your program did not run со, check fODNEBBOESE 
&TYPE Do you wish to XEDIT the program file? (Y) 
&READ VAR &RESP3 

XIE SRESP3 NER S P 2 

&COMMAND XEDIT &FN FORTRAN A | 

“ТҮРЕ Do you wish 5о run the program again? (Y) 

&READ VAR &RESP4 

&IF &RESP4 EQ Y &GOTO -H 

ЧЕЛІТ 2 

&IF &OUT EQ Y &GOTO -REDO | 

&TYPE Your output is τη επ EIN DIM κο 

“ТҮРЕ Do you wish to BROWSE your output? (Y) 

&READ VAR &RESP 

«IP &RESP EQ Y &COMMAND BROWSE &FN OUTPUT А 

«ТҮРЕ Print vour$outbput$e lle 

&READ VAR &RÉSP?7 

A EQ 7 &COMMAND PRINT «ЕМ ООТРОТ А 

= О 

&TYPE Do you wish to XEDIT the program file? (Y/N) 
&READ VAR &RESP5 

&IF &RESP5 EQ Y XEDIT &FN FORTRAN A | 

&TYPE Do you wish to run the program again? (Y) 

&READ VAR &RESP6 

«ТЕ &RESPÓ EO Y. AND ӨБЕЗР5 ЕО 775 59-ы 

&IF &RESPÉ EQ Y &GOTO -RUN 
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APPENDIX F 
UNIT TYPE AND FORMATION 


Table 18 shows the integer value and the unit type it represents. 


TABLE 18 
Meee OF UNIT 


Integer Value Type of Unit 
1 
2 
3 


Tracked Vehicle 
Wheeled Vehicle 
Dismounted Troops 





Table 19 shows the integer value and the unit formation it represents. 


TABLE 19 
UNIT FORMATION 


Integer Value Unit formation 


Deplo 


( ар Іт column ) 


Undeployment 
( μα“ ma ) 
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APPENDIX G 
RATING CONE INDEX 


Figure G.1 shows an example of how the rating cone index can be used [Ref. 4: p. 
20]. 


MAXIMUM MAXIMUM 
SLOPE SLOPE . 
NEGOTIABLE NEGOTIABLE 
(PER CENT) (PER CENT) 


10 
ы IRACKEO VEHICLES WITH Ναῖ--... € I — mM 
suu. D, LT T O 
ια а 


LESS THAN I m O 


ete 


Criteria for self-propelled vehicles on fine-grained 
soils and sand with fines, poorly drained. 


NOTE. % Slope Negotiable" 
increases little above 50% 
when soil strength increases 
from +50 to +100 


RATING CONE INDEX-POINTS ABOVE VEHICLE CONE INDEX 





Figure G.1 Rating Cone Index. 
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